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ABSTRACT \ 


The eee hy of a Goemoutational ly simple yet 
cryptologically robust rule for generating pseudorandom 
bitstreams cannot be overstated. In most applications we 
strive to detect and avoid chaotic behavior; here we embrace 
a Mee ed cilar chaotic discrete dynamical system to exploit its 
use as a driver for a pseudorandom number generator. The map 
from the Hénon attractor to the binary domain {0,1} proposed 
by Forré/Heyman has been tested cryptologically and 
Seaeastically with mixed results. In this thesis we 
mathematically evaluate this symbolic dynamics scheme and 
investigate more rigorously its utility as a pseudorandom 
number generator. Specifically, we demonstrate how the 
property of being one-to-one holds, but that the property of 


being onto does not. 
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I. INTRODUCTION 


A. HENON’S DISCRETE DYNAMICAL SYSTEM 

Michel Hénon (1931- ) waS a French mathematician and 
astronomer who was interested in modeling astronomical 
behavior. Henon learned how the meteorologist-mathematician 
Edward N. Lorenz (1917- ) had developed a three-dimensional 
attractor to model the complicated dynamics of thermal 
convection. An attractor is characterized by the following 


Serinlecion. 


Pee UNL TiLON: Let V be a subset of R” and F: V7 R", where 


[ewe 3s. Let A be a subset of Ve Then A is an attractor of 


1. A ioe meses invariant subset of V under F. The set 

fern R” 1S invariant under-F: VR" if F(A) ¢ A. 
11. There is a neighborhood U of A such that if v isin YU, 
F°(v) ~ A as noe ,where F" denotes the n-fold composition of 


feeeer. | 7p.147,200j. 


The Lorenz attractor 1s generated by a system of three 
differential equations and relies upon integration which does 
not lend itself to timely and accurate computer calculations. 
[Ref.2] Hénon sought to build a simple discrete dynamical 


system that would retain the qualitative properties of the 


continuous Lorenz system. His goal was achieved through study 
of the Poincaré section of the Lorenz equations. [Ref.3] 
Thus, the Hénon attractor 1s, to some extent, a two- 
dimensional version of the Lorenz attractor. The Hénon map is 
the following quadratic recurrence which maps the Euclidean 


plane into itself. [Ref.7] 


H : R* +> R? 


‘ 
T 
oy 

x 
pes 
Ny 


Each point in R* has a unique image and preimage under H. The 
Same is true of the inverse of the Hénon map. The constants 
a and b are real parameters. We can write the map as a 
EransLormat lOnedn ene Lo Wlovwinc foun: 

His y) =) (0 =] ax Ba eee 


Understanding the dynamics of the map 1s easier if we 


decompose it into the following independent transformations: 


H, (x,y) = (x,y + 1 - ax?) 
Bae) =— (v4) 
Hey ee 
such that 
Hx, V) sa wee ee 


Respectively the transformations cause bending, contraction 
gen) © |<1), and a reflection. [Ref.5:p.662] This stretching 
and folding are analogous to the kneading or mixing common to 
chaotic systems. 

The parameters (a,b) control whether the system forms an 
interesting attractor with apparent aperiodicity oor an 
uninteresting attractor of low period. Reportedly, Henon 
chose his classical parameters (a,b) = (1.4,0.3) because they 
gave the "nicest and strangest" picture. Figure 1 shows the 
computer representation of the Hénon attractor FOr 


more. —=(0,0) using 2000 points. 





X-axis 


Figure 1 Hénon attractor 


In keeping with our definition, the computer representation 1s 


generally unaltered by the choice of (x,,y,). This 1s true as 


long as the values tend to the attractor and the first) eee 
100 iterates are ignored. The first points generated by the 
map may not be close to the attractor, but the generated 


points get closer and closer to the invariant set. 


DEFINITION: Let x belong to the domain of a function feu 


of period n if f"(x) = x, and is of prime period (na 
x,f(x),£*(x),...,£ (x) are distinct... Ti x has. neriede ae 
the orbit {x,fi(x),f*(x),...,£7'(*%)} 1s a pericdic Borman 
[Ref.l:p.21]). If n =1, then x 1s said to be a fixed poanmes 


The merit of Hénon’s choice of (a,b) parameters can be 
found by observation of the respective bifurcation diagrams. 
A bifurcation diagram gives us information regarding the 
periodic and nonperiodic orbits of amap as a parameter 1s 
altered. In a bifurcation diagram, the vertical axis 
represents iterates of the variable x as a parameter 1s varied 
on the horizontal axis. [Ref.1:p.52] If we take a vertical 
slice at a particular parameter we get an idea of the 
attractor’s structure for that parameter. In Figure 2 the b- 
parameter 1s held fixed at 0.3 and the a-parameter is varied 


Lrom.0..2 ene toe 


n 
én 
= 
cer 
Send 
© 
x 


0.6 Oye" 
a-parameter 





Figure 2 Bifurcation Diagram: b-parameter fixed 


Notice that from about a = 0.1toa= 0.4 we have a point 
BieemaelOr, that 1s, the entire attractor consists of a single 
Mees point. From about a = 0.4 to a = 0.9 we have a period- 


two attractor; that is, the attractor 1s made up of two 


Series . These period-doubling bifurcations persist as the 
parameter a 1S increased, until a critical value of 
approximately a = 1.1 is reached, where we enter the chaotic 


regime. Notice that Hénon’s classical a-parameter 1.4 lies in 
such an area of wildly complicated dynamics, sometimes called 
a chaotic band [(Ref.6]. 

In Figure 3 the a-parameter is held fixed at 1.4 and the 


b-parameter is varied from 0.0 to 0.32. Again, Hénon’s b- 


value 0.3 not Surprisingly lies in an area Of Chace. 


aACeaiveney . 


n 
~ 
_ 
om 
i) 
oe) 
yd 


‘Seles: 8% 
b-parameter 





Figure 3 Bifurcation Diagram: a-parameter fixed 


B. NATURE AND STRUCTURE OF THE HENON ATTRACTOR 

Because the Hénon transformation is quadratic, it 1S 
possible that even relatively small (x, y,) values can produce 
an iterated sequence, or orbit, that escapes to infinity. 
[Ref.4:p.664] The Hénon attractor has an associated basin of 
attraction of points (found numerically) which in forward 
iteration are drawn to the attractor. A point in the basin of 
attraction thus tends to the invariant set. The Henon 


attractor also has an associated area Known as the trapping 


region. The trapping region is a quadrilateral within the 


basin of attraction defined by the vertices: 


Pe Na See 2) B eee 0213 3)) 


C) 
i 


(ie eeae ae 4) D 


eee. 0). 


If (x.,y,) 1s chosen from this region, all subsequent iterates 
will remain in the region [Ref.5:p.664}]; hence the trapping 


Meagvon 1S invariant. Let 


fo nenen adtberaccor 
Bee Soc imvOr “Aaberaction 


ie (bap pangs red 20n.. 


Then, we have that 

1, ac oe 
foes relationship is shown explicitly (excluding the trapping 
region) in Figure 4; the shaded region is the basin of 


Serraction. 





Figure 4 Basin of attraction (shaded 
region) and Henonmateracten. sus 
figure is from reference 5S, p.605. 


Therefore, 
1. if (X,,Y,) © UR, then (x-.v.) ¢ IR Om senor 
11. 1f (x,,y,.) € B, but (x ,.y.) € TR, then (<9 7.) may gee 
belong to TR for n < N. However, from the description 
of the basin of attraction we know that for all n 


sufficiently large, si yewe Th. 


The b-parameter, 0.3, is in the dissipative term of the 
relation. The Jacobian of the Hénon transformation is (-b), 
and hence the map is dissipative when |b| < 1. (Ref.1:p.168]. 


Therefore, should we iterate the entire trapping region under 


MiemesaistOrmation for b = 0.3, this area would shrink and in 
Memlimit tend to zeros Lieeomem cee ractor thus has. zero 


area. Shown symbolically, 


Pee ee ie. eee ee) 


ou 


Mi@eract, any point belonging to the trapping region can be 
iterated to construct a dense orbit on H,, except the periodic 
points which are countable and hence have Lebesgue measure 
zero. The visual result, in the limit, is the same for any 
initial condition (less a finite number of initial iterations 
memeassutre close proximity to the attractor). Two different 
points will generate different sequences but the limit set, 
under the above qualification, will appear to be the same 
visually. Furthermore, in general, there 1s no correlation 
between the trajectories of two different initial points, 
except possibly inthe first relatively few iterations should 
the initial points be neighbors. (We disregard here the 
instance where, say, trajectory 1 has common points with 
Mmeamectory 2 but for different iterates. In this case the 
orbits are shifted versions of one another.) This lack of 
correlation is related to sensitive dependence on initial 


menaitions. 


The Héenon equations (1.1) and (1.2) (with the classmeam 

parameters (a,b)=(1.4,0-.3)) May be wrieotmeerre variable as: 
Xp Clad 0 eee (lee 

There are two fixed point solutions to the recurrence, where 

ea = ee 


n+} n 


This occurs when (x,y) = (.6314...,.1894...) (fP2) andes 
= (-1.1314..., -.3393...) (FP1). The fixed points, txapeeeme 


region and Hénon attractor appear in Figure 5. 


oo wongitudinal Axis 


Transverse Ax13 





Figure 5 Fixed points of Henon attractor, trapping 
region and eigenvectors associated with FP2 


The left fixed point FPl is repelling and 1S noticeably not on 
the attractor or even within the trapping region. The 


remaining fixed point FP2 is a saddle point and belongs to the 


Lo 


aeeractor . The Jacobian matrix at FP2 has associated 
eigenvalues A, and A,, with slopes p, and p, for the respective 
eigenvectors: 

ye 01559... pion 37 eS. 

Jay, ES ee) OA ie ee a 
femme 5 also shows us the strong local effect of the 
Srgenvalues at FP2. Longitudinal spreading is due to A,. 
Since |A,| is greater than 1, points are forced outward from 
FP2 along the longitudinal axis forming a roughly broken line, 
at least locally. The eigenvalue A, causes a transverse 
@emeraction of points. Figure 6 shows a detail of this 
Structure. The dense transverse structure is revealed through 


eee) Magnification in Figure 7. 





Figure 6 Magnification of area around fixed point 


la 





Figure 7 Detail of dense transverse structure around 
fixed point 


As magnification 1S increased and more points are plotted the 
attractor begins to take on the appearance of distinct bands 
with a cross-section that has the nature of a one-dimensional 
Cantor set -- a totally disconnected set of points with 
patterns that persist on infinitely small scales. This is the 
self-similar nature of the attractor. Thus, at least ona 
level local to FP2 we observe a fractal quality that 1s common 


in chaoticta@eractors (Ret > -p 67 ols 


C. CHAOS 
Because the field of chaos is in its relative infancy not 


all terms are universal. The most mathematically rigorous and 


eZ 


widely accepted definition of chaos for discrete dynamical 


systems is the following: 


Pee iTtlONe Let S be a set. A mapping f: S—3 S is said to be 
eractic on S if 
1. £ has sensitive dependence on initial conditions, 


fm ft 15 topologically transitive, 


LJ 


the set of periodic points is dense in S. [Ref.8:p.50] 


We will now elaborate on each of these concepts ina bit more 
detail starting with sensitive dependence on initial 


Sendai c1ons. 


DEFINITION: A mapping f: S 4 S has sensitive dependence on 
initial conditions if at every xe S, there exists an € > 0 
Slelethnat for each 6 > 0, there is a y in S and a positive 


integer n such that 


lx-y |< 6 and Fee re neers var | sno Lee eel. 64 


Based on empirical evidence, the Hénon attractor is widely 
believed to possess sensitive dependence on initial 
conditions. [Ref.7] The stretch-and-fold phenomenon mentioned 
previously ensures that close points will not be neighbors 


meee 4a f£inite number of iterations. Figure 8 shows a plot of 


3 


the Euclidean distance between the orbits of two “2niseeeem 


COnaditLons: 
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Figure 8 Sensitive dependence on initial conditions 


The two points are the origin (0,0) and a point (0,€’) where 
€’ is just within machine epsilon, that is €’ < €. Machine 
epsilon in this case is approximately 2.22 x 10°°°. Within 
seventy iterations the orbits are completely different. 

Of course, since both initial conditions belong to the 
trapping region the maximum difference between the orbits 15s 
bounded by the extremes of the trapping region. Although this 
1s only the first of three criteria for chaos some accept the 
dynamics on the Hénon attractor as chaotic based on sensitive 


dependence alone. [Ref.4:p.671] Sensitive dependence brings 


14 


out an interesting paradox. Globally we have an attractor 
Peele reSspece to the basin of attraction; however, locally for 
the points of the invariant set, the orbit jumps around the 
attractor chaotically, seemingly repelled by each point. 
Memelogical transitivity, also called mixing, ensures that 
the set on which the dynamical system is defined cannot be 
Seeeoed into two disjoint open sets which do not “interact" 


(Or more precisely, are invariant) under the system. 


DEFINITION: A mapping f: S—~S is topologically transitive if 
for any pair of open sets U, V, which are subsets of S, there 


peeeeeoek > 0 such that £*{U) AV #0. [Ref.8:p.49] 


Essentially, for a system to be topologically transitive 
memmust be true that for any two arbitrarily small regions U 
and V (where U and V are not points) we can find an initial 
point x, 1n U whose orbit will enter the other region at some 
Meerarion. (Ref.5:p.554] 

For a system to possess dense periodic points means that 
Mmaimeitaroitrarily small neighborhood of any point belonging to 
the attractor there must be a periodic point. This condition 
Suggests paradoxically that one of the qualities for "chaos" 
1s a high degree of "order". 

Although there exists vast numerical evidence to support 
the claim that Hénon’s dynamical system possesses all three of 


Mmiyese Criteria for chaos, none of the criteria have been 


i) 


rigorously proven. [Ref.5:p.671] [Ret.19:p.152] To kéepiaa 
in perspective, it must be acknowledged that our definition is 
more stringent than proving the existence of a dense orbit and 
a positive Lyapunov exponent, a commonly accepted alternative 


Gefingeaon. hema 


D. THE SHADOWING LEMMA - TRUSTING THE COMPUTER 

The study of theoretical chaotic dynamical systems would 
be extremely handicapped if we could not trust our computer 
representations of these systems. Because a computer 15 
limited in terms of its accuracy by rounding errer eee 
appropriate for us to ask whether our computer calculation 1s 
an approximation of some true orbit of the system. [Ref.9] 
The Shadowing Lemma shows that we can trust the COMmpUres 


image. 


DEFINITION 1: The map f£:R"°"~ R is a class C® map if it 1S p 
times continuously differentiable with respect to x,, x 77 eee 
for some p, 1S pse. A mapping £:R"—- R" is class CY’ ltveaen 


component f, of f, (where i=1,...,m) 185 p times coOntinvgelia, 


differentiable [Ref.10:p.1-2]. 


DEFINITION 2: A mapping f:R"°—> R" is a diffeomorphism if it is 
one-to-one, onto and both £ and £+* are continivemem: 
differentiable maps; f is a C®’-diffeomorphism if both £ and f° 


ane C* maps [Reiakt >. 2e 
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DEFINITION 3: Let f: R* 4 R* be a diffeomorphism. A set A is 
said to be a hyperbolic set under f if: 

mares oll pornts pp belonging to A. there is a Ser of 
lines E*(p) and E*(p) in the tangent plane at p which are 
preserved by the Jacobian of f at p. (E*(p) is called the 


M@eeample line and E’(p) is called the unstable line). 


meee Oo) momG FE ip) vary Continuously with p. 
3. There is a constant A>1 such that |J(p)(v)|2Alv| for 
peeve (p) and |J“*(p) (v)|2A|v| for all veES(p). [Ref.8:p.236- 


Ze! | 


Suppose f satisfies the three previous definitions for some p 
Wath an associated hyperbolic invariant set A. The true 
(forward) orbit of a point z, under the map f£ 1S given by the 
Seemence {Zz}, such that z, = £f'(z,). [Ref.10:p.351] Due to 
finite precision the exact calculation of the sequence 
{f£*(z,)}-" 1S impossible. Because we are unable to obey the 
memeerton Zz = f*(z,) in order to find an exact orbit, let us 
instead define a pseudo-orbit or more specifically €-pseudo- 


ee) t.. 


DEFINITION 4: An €-pseudo-orbit is a sequence of points, 


1Yy,}- , such that y. € A, and 


Oy ate ye 


es 


where dis a metric on R". ([Ref.10:p.352] The €-pseude-Graae 


1S PlcEured at EMenrop of ferguce 9. 





Figure 9 €-Pseudo-orbit and 6-Shadow 


Notice that we adjust each iterate by dropping f(y,) and then 
choosing y,,, to be within ¢ of f(y,): We can see thatiieme 
term pseudo-orbit, therefore, 1S appropriate. One final 
ope alighit & 1edo 1s necessary for us to understand the Shadowing 


Lemma. 
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Magee TTON S: Suppose {xX,}, 1S an actual orbit, that is, ~xeA 
e@eex, = f*(x,); {x,}> 1s a O-sWadow of the pseudo-orbit {y }-* 
if 

iG ures 0 
The 6-shadow of the pseudo-orbit is shown at the bottom of 
Smoure 9. The Shadowing Lemma assures us that within an 
arbitrarily small distance, 6, of an €-pseudo-orbit there is 


ammpemact orbit. Stated succinctly: 


Shadowing Lemma: If Definition 2 holds for the map f and Ais 
a hyperbolic invariant set on R", then for every 6 > 0 there 
1S an € > 0 such that every €-pseudo-orbit in A is 6- 
S@eeowea by the actual orbit of some point x eé€ A. 


Weert. 10:p.352)] 


Although the Shadowing Lemma is very powerful, 1t does not 
explicitly speak to us about the computed orbit, call it w, 
Bor which we hope to find to an exact orbit nearby. 

®@ur Computed orbit, w,, 18 not the true orbit {f (y-)}, we 
are looking for; however, let us show more clearly why it 1s 
an acceptable approximation of a true orbit in the system. 
(Ref.9:p.251] The sequence of computed values, w., from the 


Henon map for some initial condition y. is actually the 


floating point representation of f(w..,) or, 


IS, 


Wa. (= 6 


W, = eieeio ee 

w, = E1(f£(w,)) 

Wap = ere ee oe 
Recall from our definition that f(ly,) 1s within € Of yee 
order for us to show that a computed orbit, w;, can be 


arbitrarily close to an exact orbit, x,, which 6-shadows the 
orbit, y,, we must first make an important assumption. We 
assume that at each iteration the error between the computed 
value of fly,) and the exact value of fly,) 1s bounded. That 
SG 

d(w , £ (ysee eee such theipmicumerc, 
where €, is some function of the machine precision which 
bounds the ‘eryvwer" for al ieeeecrarions. Figure 10 shows a 
representative configuration of the orbits involved. By the 


Shadowing Lemma d(y,,x,) < 6. 


oN 





Figure 10 Pseudo-orbit, exact orbit in 6-shadow, and 
computed orbit 


Using the triangle inequality it follows that 


CMe eel (oe ly) ) er Kw, Pty.) 
mere , 
Meoreety ei S Gi xk ay e+ ai tty.) )-. 
Therefore, 
Giweree = “dix yi) “= adty.,f(y.3)) = diwatiy.)) 
or 
clique < iS a) + E + E.. 
If we let 6’ = 6 + € + €, then it follows that, for every 0’>€, 


we can calculate a computed orbit, w,, that is 0’-shadowed by 
an exact orbit, x,. In this way we see that the existence of 
our computed orbit in the 6’-shadow of an exact orbit actually 


depends on the precision of the machine in use. 


Zt 


According to a theorem in a recent paper (Ret Ji5)5 
specific criteria have been developed under which we are 
assured that in the shadow of a computed orbit there will 
always exist an exact orbit within a calculable tolerance. 
The authors accepted that the Hénon map and the Héenon 
attractor satisfied Definitions 1 and 2 previously mentioned. 
The accuracy of the numerically computed orbits which the 
theorem measures applies to non-hyperbolic sets as well as 
hyperbolic sets so compliance to Definition 3 may be relaxed 
[Ref.17]. Computations applying the theorem in reference 15 
were performed in Microsoft Quickbasic using a double 
precision IBM system computer. The value of € (the distance 
between y, and f(y,_,) which defines the pseudo-orbit) was held 
to approximately 2°*°. The authors applied the theorem a@acuus 
was discovered that even after 372,000 iterations there is a 
true orbit which differs by at most 2° (approximem see 
1.8626x10°) from the computed orbit generated by the Hénon 
map. 

Based on the dense structure and infinite detail of the 
Hénon attractor it is intuitively reasonable that in the fooer= 
Print of the computed orbit there exists an exact Wememan 
Thus, the Shadowing Lemma assures us that the statistical 


evidence measured under computer analysis is significant. 
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II. PSEUDORANDOMNESS FROM CHAOS 


A. LINKING CHAOS TO THE BINARY DOMAIN 

Symbolic dynamics 1S a technique that can be used to 
relate the dynamics of a particular system on a metric space 
to another system on symbol space. Symbolic dynamics allows 
us to analyze a system by studying its effect on symbol space. 
The principal idea is that if two systems are topologically 


conjugate, then their dynamics are equivalent. 


DEFINITION: Suppose we have two maps: f£: U7 Uandg: VV. 
The functions f and g are topologically conjugate if there 
exists a homeomorphism h:U-7V such that hofz=gqoh. The 
function h is a homeomorphism if: 
1. his one-to-one 
MEEINETION: Suppose £f aS a function from U to V; £ is 
one-to-one if and only if for all elements u. and u.-éeU, 
Mme) =r) dlimolres Uy = u-. 
meet) 1s onto 
DERG TION: Stieeose £f is a titimction from U to V; £ is 


onto if and only if for any element v in V there 


Ssetcecwan element’ u in U such that v = £(u). 
mer. ff 1s continuous 
j/mem hh’: is continuous. 


Zo 


Criterion iv is actually implied by criteria 1, 11, ance 
Furthermore, if £ is topologically conjugate to g €tHemuume 
follows that ¢@ Ts" topedloqg testy cengigaeon ee mas 

Finding a homeomorphism that forms a topological conjugacy 
between two maps can be extremely difficult. However, this 
does not diminish the power of the relation. A topolo@qmaaa 
conjugacy, as we mentioned, relates the dynamics of two 
systems completely. [Ref.20:p.27] For example, i1f the 
function £, from our definition, has a period-two cycle wigmaae 
then {h(p),h(q)} 1s.a period-two cycle for the Eunct 1 Ommgeaen 
this way, all orbits for the function £ have correspomgamae 
orbits for the function g. Furthermore, if £ has a densemeee 
of periodic points in U, then the same is true for g in V. 
[Ref .5:p.571] Following this reasoning, if it can bew@eitern 
that the dynamics of map f exhibit the three characteristics 
of chaos previously mentioned, and there exists a topological 
conjugacy between f and g, then map g will also be chaotic. 
[Re fsa pec, 

Can we harness this chaotic energy that exists in one 
system and through an effective symbolic dynamics transfer 1t 
to another under a homeomorphism to achieve a useful result? 
This 1s precisely the issue raised in Réjane Forré’s treatise 
of November 1990. [Ref.12] Forrée hoped to apply a presumably 
chaotic discrete dynamical system to the field @iex 
Cryptography. Specifically, she attempted to devise a scheme 


to generate nearly random sequences of zeros and ones which 
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could be used for coding purposes. (Any numbers generated by 
a deterministic rule cannot be truly random and are, 
therefore, termed pseudorandom.) Forré’s proposed symbolic 
dynamics relates the apparent chaotic dynamics of an orbit on 


the Henon map to binary codespace by 


h:H, 2 d.. 


Here, A represents the attractor associated with the mapping 
mumps = Gefinition Chapter I. section A. p.1); in particular, 
f£ is the Hénon map. The set »,, also called symbol space or 
codespace, represents the collection of all infinite sequences 
of zeros and ones. Any chaotic behavior exhibited in the 
Hénon map transferred to }, would be observed as a 
pseudorandom stream of zeros and ones in 2... Forré’s symbolic 
dynamics is based solely on the horizontal component, x,, of 
the iterates of the Héenon map. The elements of such a binary 


N 


sequence {S,},.," 1n 2%, are defined as follows: 


vie ie ene Se =. Os 


I 
ke 


iG oe eee LEI yo) 


Mer a string of length N where x,.. represents a dynamic 
median. (The dynamic median ensures that on the average the 


Erajgectory will fall on each side of the median for half the 


ZOD 


iterates; this median has been more carefully calculated as 
Xyon=- 4098 iteewetereuies maken In order to clarityaieae 
computational process involved in Forré’s symbolic dynamics 
the following example is furnished. Let us arbitrarily choose 
a particular (x,,y,) value from the trapping Yeqivenee 
.(x%,,y.) = (1,0). Since x, > X. the binary sequence Gilemem: 
corresponding to (x,,y,) 1S S, = }. Using equations (335i 


(3%2) below we calculate Vy was. Bemtavce 


X,.1 = 1 - 1.4 x,’ + y, (ee 
inet BPG (3 an 
a eer ry ores — 2 


Because xX, S Xyp the next binary sequence element, 
corresponding to. y,) 7s Siem In this way we can 


calculate the entire forward sequence of binary elements: 


S = (So; 04 poo oO. ee = (ele Ourde Occ.) 2 


Although Forré’s calculated median was inaccurate (she 
used a dynamic median x,,, = -39912) her results showea@iaa 
under the cryptographic properties of linear complexity and 
Jump complexity the bitstreams were wholly indiscernible from 


truly random sequences. However, she conciuded that a third 
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mipepercy, the n-tuple distribution, was inconsistent with 
truly random sequences and that this rendered the unaltered 


scheme unsuitable as a pseudorandom number generator. 


Be RECENT STUDY OF FORRE’S SYMBOLIC DYNAMICS 

Ferre Grew attention to the poor n-tuple characteristics 
of her unaltered symbolic dynamics. In amore recent study 
feet os} Heyman refined the calculation of the median, x,,,, 
and rigorously investigated the n-tuple property, as well as 
three other cryptographic properties, to evaluate the 
pseudorandom number generator. It was concluded that the n- 
Spee Or runs property was a minor detractor and that the 
proposed use of the scheme as a cryptographic pseudorandom 
number generator was sound and effective. The runs property 
was mentioned by both Forré and Heyman, but the authors’ 
Semelusions with regard to the significance of the property 
were markedly different. Under these mixed results, the runs 
anomaly demands further investigation. 

Due to the disparate conclusions drawn by Forré and 
Heyman, it 1S a natural next step to evaluate the symbolic 
dynamics under more mathematically rigorous criteria in order 
to determine more objectively and conclusively whether the 
scheme is an effective pseudorandom number generator. Because 
the chaotic behavior on the attractor and certainly the 
symbolic dynamics scheme itself are mathematical concepts, we 


evaluate the pseudorandom number generator mathematically. 


ou 


This 1s certainly necessary since the previous papers were not 
analytical and gave only a cursory mention of the theory of 
chaos. Primarily, our task is to collect evidence to prove or 


disprove that this scheme, 


heh, = 


gives a homeomorphism. Our secondary objective is to explain 
completely the runs anomaly in mathematical terms. Thuse 
based on the definition of a homeomorphism we will provide 
experimental evidence to support the one-to-one property. 
Furthermore, we will provide analytic proof that the presumed 
homeomorphism 1s not onto for the proposed pseudorandom number 


generator. 
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III. EXPERIMENTAL EVIDENCE FOR ONE-TO-ONE 


A. A COMPUTER PARADIGM 

Mine Henon attractor A = H, 1S an infinite set of points 
feeecuced by the Henon map. The initial condition of the Hénon 
recurrence dictates which orbit on H, we produce. Let A _,A- be 
two infinitely long sequences of points (orbits) based on two 


different initial conditions (s,,t,) and (l,,m.) respectively, 


ieee {S_/t,), (1,,m) are from H,. The sequence A. will be 
Memmmed as {(s,,t,), (S,,t.), (S:,t.),.--}4 The sequence A. will 
eee ene), (1.-,m,), (1,,m,),"..}.~ Although the two-dimensional 


pmoes Of A. and A, are indistinguishable, the two sequences of 
Pemmes A, A, are not the same. Let h(A,),h(A.) belong to Jy, 
the space of all infinite sequences of zeros and ones. eng 
the symbolic dynamics h:H, ~ 2, to be one-to-one, it must be 
Mem@emechat if h{A,) = h{A,) then A, = A. That is, if two 
binary sequences belonging to »}. are identical then they must 
Beetiap oped from the same initial point (x,y) in H,. 

We cannot hope to find with any precision the points 
Mmewenging to H,. Furthermore, by increasing the number of 
possible initial conditions or keys from which we can generate 
@ietmes Close to the attractor we enhance the cryptographic 


qualities of the pseudorandom number generator. That 1S, we 


reduce substantially the possibility of the key being found 


29 


and the code being broken. For these reasons we modify our 
symbolic dynamics scheme to h:TR—~%,. That is, we increase 
the number of possible keys to all those possible points from 
the trapping region. Under the Hénon map the trapping region 
is also invariant; therefore, h:TR~», is a symbolic dynamics 
analogous to h:H, > d». 

If two binary sequences are not the same then they must 
have originated from different initial (x)y) valwess The 
total count of computer representable numbers in the trapping 
region, although extremely large, is finite. Thus the task of 
proving that the homeomorphism 1s one-to-one on an uncountably 
infinite number of points in the plane is avoided and 
unnecessary. Certainly, no irrational numbers are computer 
representable. Since all pseudorandom number generators are 
implemented in a computer environment we qualify our goal to 
showing that the presumed homeomorphism 1s computer one-to-one 
on TR. That 1s, we attempt to show that the homeomorphism 1s 
one-to-one with respect to the computer representable points 
in the quadrilateral. 

We propose then a computer paradigm that models the finite 
number of computer representable numbers. The model will also 
give us insight into the complicated dynamics of the map. 
Consider the trapping region or quadrilateral as a grid 
consisting of a large but finite number o& points. Figqguiemam 


shows the quadrilateral and the median x,,, = .4098. 
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Figure ll Trapping region and dynamic median 


The quadrilateral is divided into two sides by the median. 
An initial point chosen from the left quadrilateral will 
correspond to a sequence under h:TR—~ 3%, that will begin with 
the binary digit zero (0). Similarly, an initial point from 


Bme Fr2ght quadrilateral will correspond to a sequence tEnat 


PoeeeeSegin with the binary digit one (1). Let us focus om the 
right quadrilateral. (Our argument will extend to the left 
quadrilateral.) 


We model the computer points in the right quadrilateral 
region by first enclosing it in a rectangle of minimum area. 
Byechoosing a particular spacing we can fill the rectangle 
with equally spaced points to produce a Cartesian coordinate 


system. If we disregard the points in the rectangle but 


By 


outside the trapping region we have a grid which can be used 
to model the large but finite number of computer representable 
Peles: If we can show conclusively that as our Spaegmg 
becomes smaller and smaller, the binary Sequences of mmee= 
points are all different from that of a point chosen at random 
not belonging to the grid points, then we should be convinced 
that h is one-to-one. Fundamentally, we first choose a fine 
grid spacing for the right quadrilateral that correspendemiae 
a large number of points inthe plane. As we iterate all the 
grid points and our random point through the Henon map weseuim 
out those points that at each iterate do not give the same 
binary element from {0,1} as the random point. We must 
convince ourselves that there exists some iterate for which 
all the binary sequences corresponding to the grid points 
differ from the binary sequence of the random point. Thus, we 
focus on a succession of subsets of our original Grid @emmee 
at each iteration. In this systematic way with computer 


assistance we hope to reveal the answer to our question. 


Be PUTTING THE MODEL TO WORK 

Let us begin with an example of our model in practice. This 
example was chosen because it allows the reader a clear, 
typical depiction of the procedure in a small number of 
iterations. Various MATLAB programs for this process on Bama 
sides of the quadrilateral can be found in Appendix A under 


the names GRDCOMP1.M through GRDCOMP4.M. As shown in Figure 


BZ 


izpene model lets us Ei1ll the right side Of the quadrilateral 
with equally spaced points determined by a specified value for 


the spacing. 


Grid Field and Pointe Chosen on Attractor 





Figure 12 Typical run right quad model: before iteration 


Here we use a relatively coarse spacing of 0.06. The point 
(W,Z) represented by a cross within the field is chosen at 
random but should be different from the grid points. Figures 
13 through 21 represent subsequent iterates of the Hénon map 
where only those points from the grid which give the same 
binary sequence as the randomly chosen point (W,Z) are 
preserved. The median, xX, = .4098, 1s the dynamic median, as 
previously mentioned; therefore, it 1S not Surprising that we 


lose approximately half the points at each iteration. 
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Movement of Pounss in Grid with Same Binary Sequence 





Figure 13 Typical run right quad model: Iterate l 


Movement of Points in Grid with Same Binary Sequence 





Figure 14 Typical run right quad model: Iterate 2 
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Movement of Poinss in Gnd with Same Binary Sequence 





Figure 15 Typical run right quad model: Iterate 3 


Movement of Points in Grid with Same Binary Sequence 





Figure 16 Typical run right quad model: Iterate 4 


Les 


Movement of Points in Grid with Sarne Binary Sequence 





Figure 17 Typical run right quad model: Iterate 5 


Movement of Points in Grid with Same Binary Sequence 





Figure 18 Typical run right quad model: Iterate 6 
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Movement of Pounes in Gad with Same Binary Sequence 
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Figure 19 Typical run right quad model: Iterate 7 


Movement of Poinss in Gad with Same Binary Sequence 





Figure 20 Typical run right quad model: Iterate 8 


oy 


Movement of Points in Gnd with Same Binary Sequence 





Figure 21 Typical run right quad model: Iterate 9 


By iterate 10 (not shown) no grid point is on the same side of 
Xyen @S (W,Z). Figure 22 gives us an idea of how quickly the 
Bolnt S Ehat follow the orbit of (W,Z) Giminish in number. The 
curve (1/2)" 1s plotted as a comparison. TABLE 1 shows the 
count of points for iterates 1 through 9 that give the same 


binary sequence as (W,Z). 


TABLE 1 NUMBER OF POINTS THAT GIVE THE SAME BINARY SEQUENCE 


a 


NUMBER OF POINTS THAT GIVE THE SAME BINARY SEQUENCE 
mare 0 err a ee 
NO. PTS 98 65 13 g 4 | 4 4 4 4 1 
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Decrease in No. of Grid Poincs with Same Binary Sequence and (4)" vs. n 
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Iteration 
Figure 22 Grid points with same sequence after n 
1terations 





Presumably based on the map’s sensitive dependence on initial 
conditions there occurs a natural spreading effect, most 
noticeable in the left half of the quadrilateral, that forces 
the previously close points away from each other. A typical 
example of a run of this procedure using the left 
quadrilateral can be found in Appendix B. 

Does there always exist an iterate N where the sequence 
for (W,Z) diverges from all sequences corresponding to even an 
extremely fine grid? Certainly this example only gives a. 
Mase’ OL the capability of the computer to fill the right or 
mee Guadrilateral with points. Based on the limitation of 
computer memory it is not possible to use a Cartesian spacing 


small enough to model even a small fraction of the huge number 


oe, 


of computer representable points. This is not our goal. We 
concern ourselves, instead, with.the trend in the iterate N 
where all the grid points diverge from the iterate of (W,Z). 

Let N be the iterate at which all the sequences 
corresponding to the grid points diverge from the sequence 
corresponding to (W,Z). Using our last example we plot the 
grid point that lasted until iterate N (call it (177) 
(W,Z) for the 9 iterations. The trajectory of the two points 
1s depicted in Figures 23 through 32. Point (W,Z), agaaiuuee. 


depicted es the cross (+) and the point (u,v) by an (x). 


W,Z and Grid Point with Same Binary Sequence 





Figure 23 (W,Z) and (u,v) prior to iteration 
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Movement of Pt Chosen on HA and Marching Seqoence Pt unal Diverge 





Figure 24 (W,Z) and (u,v): Iterate l 


Movement of Pt Chosen on HA and Matching Sequence Pt unal Diverge 





Figure 25 (W,Z) and (u,v): Iterate 2 
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Movement of Pt Chosen on HA and Matching Sequence Pt unal Diverge 





Figure 26 (W,Z) and (u,v): Iterate 3 


Movement of Pt Chosen on HA and Matching Sequence Pt unnl Diverge 





Figure 27 (W,Z) and (u,v): Iterate 4 (Notice that the 
two points are essentially superimposed). 
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Movement of Pt Chosen on HA and Matching Sequence Pt unal Diverge 
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Figure 28 (W,Z) and (u,v): Iterate 5 (Notice that tne 
two points are essentially superimposed). 


Movement of Pt Chosen on HA and Matching Sequence Pr unal Diverge 





Figure 29 (W,Z) and (u,v): Iterate 6 (Notice that the 
two points are essentially superimposed). 
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Movement of Pt Chosen on HA and Matching Sequence Pr anol Diverge 





Figure 30 (W,Z) and (u,v): Iterate 7 (Notice that the 
two points are essentially superimposed). 


Movement of Pt Chosen on HA and Matching Sequence Pt anal Diverge 





Figure 31 (W,Z) and (u,v iterate « 
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Movement of Pt Chosen on HA and Masching Sequence Pt unal Diverge 
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Figure 32 (W,Z) and (u,v): Iterate 9 


Meer1ce the proximity of (W,Z) and (u,v). They are in no 
way superimposed. Indeed, ina Euclidean sense the two points 
are not close at all. These positions may seem acceptable but 
even when the spacing between points 1S dropped to the limits 
meee ompuLlcer memory the proximity of (W,Z) to (u,v) does not 
appreciably change from this typical example. It is a fact 
that throughout the trajectories of the two points there may 
be 1terations where the points are virtually superimposed as 
in this example, but as a rule (W,Z) and (u,v) are not 
Mztially neighbors. Furthermore, as in the typical run using 
the left quadrilateral in Appendix B, sometimes more than a 
Single point (u,v) gives the same binary sequence as (W,Z). 


i~emealhy one, buE up to seven points (u,v), in arbitrary 
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observations, have mirrored the bitstream of (W,Z) and these 
multiple matches likewise are not tightly packed points. This 
phenomenon lends support to our model because it 1s not the 
points closest to (W,Z) which typically give the same binary 
sequence. Typically many other points were much closer ina 
Euclidean sense to (W,Z) but they were rejected. That is, we 
cannot look only in the immediate vicinity of (W,Z) for some 
that offer “gO0Cd Misi estveam Smatenes- Certainly, “2a 
cryptographer knew that his chances for regenerating an 
identical binary sequence were much higher in the vicinity of 
(W,Z) he could use it to his advantage. Therefore, this 
quality is beneficial. 

The GRDCOMP1.M - GRDCOMP4.M models (mentioned previously) 
perform a grid comparison between the chosen point (W,Z) and 
the grid developed on a particular side of the dynamic median. 
The models GRDCOMP1.M and GRDCOMP3.M examine behavior in the 
right quadrilateral, and GRDCOMP2.M and GRDCOMP4.M model the 
left quadrilateral. The GRDCOMP1.M and GRDCOMP2.M models use 
a predetermined (W,Z) value from their respective side of the 
quadrilateral but not on the attractor (H,) to compare against 
the developed grid. The GRDCOMP3.M and GRDCOMP4.M models 
differ from GRDCOMP1.M and GRDCOMP2.M in that (W,Z) is 
iterated through the Hénon map to ensure that the point 


neighbors the attractor before the companucons feo ume 
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These computer models give us a feel for the N value to 
expect when we designate a particular spacing in the right or 
left grid field. The computer models PROOF1.M - PROOF4.M are 
identical to the GRDCOMP1.M - GRDCOMP4.M models except that 
they allow us to vary the spacing. We are able to designate 
a coarse and fine spacing range and an increment to use 
between these values. As the Cartesian spacing 1S steadily 
decreased and more points fill the particular half-field, by 
observation, N likewise seems to increase steadily. BE oistes 
behavior is pictured in Figure 33. Here, we use the PROOFI1.M 
model incrementing the spacing from .003 to .5 by .001 for a 


seman of nearly 500 entries. 
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Figure 33 Reduction in spacing vs. N (PROOF1.M model) 


4’] 


Although there are a few outliers all the N values appear to 
be clustered and to steadily increase overall as the spacing 
becomes finer and finer. This behavior was observed for the 
PROOF2.M - PROOF4.M models as well. 

Let us establish .10 as the coarse spacing upper bound 
because it is a round number and it allows us to fill either 
half-field with more than just a few points. In computer runs 
using this upper bound and steadily finer spacing (to the 
limits of computer memory) we observe the same general 
relationship between the decreasing spacing and N values. 
That 1s, as the spacing 1s steadily decreased the value of N 
correspondingly increases. 

To further test our models certainly we must not limit 
ourselves to a comparison of a single (W,Z) value over the 
range of spacing values. Because we observe the same general 
behavior in te GRDCOMP models whether we use a particular 
(W,Z) not on the attractor. or a (W,Z) from the attraceemae 
us choose an arbitrary number of (W,Z) values and repeat our 
models (PRF1.M - PRF4.M). For simplicity the arbitrary number 
of (W,Z) values (31) are taken from a line belonging to the 
particular side of the quadrilateral. Figure 34 shows a run 
of model PRF1.M (which corresponds to our GRDCOMP18M mode ame 


for multiple spacings and (W,Z) values). 
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Spacing 


Figure 34 Reduction in spacing vs. N (PRF1.M model) 


This run again is typical of all runs of the models PRF1.M - 
Peer. This run increments the spacing from .01 to .10 in 
steps of .001. Notice that there do not appear to be 31 
values fcr each spacing value. This phenomenon actually 
supports our observations because many of the points are 
superimposed on each other on the graph. Again, we observe an 
apparently steady increase in N values as the spacing becomes 
Meer Wich few significant outliers. BU GERS Bee hats ie: GUM ele lglet= ya) 
almost 3000 points. 

In order to show experimentally that the observed 
behavior exists, we select a modest fine spacing as a lower 
Beunad and a tiny increment. By using an extremely small 


Spacing we are able to collect an enormous amount of data to 
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use as statistical evidence. The limitation is not Comemee 
memory but the mainframe eee Statistical system (AGSS) 
which allowed a sample size of 45136 points. For each of the 
31 (W,Z) values 1456 separate runs were completed. This 
corresponded to an increment of .0000625 between .01 and .10. 


The results are pictured an eeaoizese so —ace 


NONLINEAR CURVE FIT 
SCATTER PLOT, SS7=45136 


Y = AxXEXPON Bxx 
WHERE: A, B = 21.092, —10.953 
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Figure 35 Nonlinear curve fit model PRF1.M 
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SCATTER PLOT, SS2=45136 
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WHERE: A, B = 24.989 
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Figure 36 Nonlinear curve fit model PRF2.M 
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Figure 37 Nonlinear curve fit model PRF3.M 
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NONLINEAR CURVE FIT 
SCATTER PLOT, SSZ= 45136 


Y = AxXEXPON Bxx 
WHERE: A, B = 20.824, -—9.8852 
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Figure 38 Nonlinear curve fit model PRF4.M 


The figures correspond to the models PRF1.M to PRF4.M. The 
nonlinear curve fit takes into account all 45136 points 
although only a representative handful are pictured. The 
curves respectively correspond to the following nonlinear 


equations: 


PRF1.M: curve l 2, 092m 
PRF2.M: curve 2 24.989 e *™* 
PRE Ge otaaome 182506 eon aa 
PREGEM: *GeuUrven4 J20ee 72 eeueas 


All four nonlinear curve fits are monotonically increas@mag 


from right to left. This result seems to suggest that ina 
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Statistically acceptable fashion, given a particular spacing 
for our grid we can predict an N value for which the binary 
sequence of (W,Z) will diverge from those of the grid field. 

How small can we make our grid spacing before our models 
mame  CCrtLainly, for any particular computer there will exist 
peeeoee. call if €, such that computer numbers within that 
tolerance will be considered the same. That is, if x andy 
are computer representable numbers and |x-y | < €, thenx=y. 

The value of € varies. Let computer epsilon be defined as 
avalue, €, below which 1 + € = 1; 386 Matlab version 3.5M 
computes this value as approximately 2.24 x 107'°. However, 
the value of € (using the same definition) such that 0 + € = 
Mupemeon the order of 10° for the same system. One can 
speculate that at some diminutive spacing S, our model will 
fail. That is, there will exist a spacing where the computer 
peeeenot be able to build the grid. The computer model will 
calculate the first point x in the rectangle which encompasses 
the quadrilateral. But, because the spacing S$ 1s smaller than 
the € tolerance that the computer needs to recognize the next 
point y, x and y are seen as the same point. Only the first 
Sema polnt will appear in the grid. However, this 1S in 
keeping with our definition of one-to-one for our presumed 
homeomorphism because 1f two sequences are wholly the same 
then they come from the same initial condition. That is, the 
sequences are generated from two initial conditions (points) 


that the computer recognizes to be the same in finite 


De 


precision although the initial conditions would differ in 
infinite precision. Thus, we have experimentally shown that 
for the symbolic dynamics i1f we have two unequal initial 
conditions then after some iteration Ni the correspemean. 
sequences should diverge. Therefore, we have experimental 
evidence to support the statement that the symbolic dynamics 


1S one-to-one. 
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IV. ONTO 


A. ONTO PROPERTY 

Whether or not the proposed symbolic dynamics has the onto 
property previously mentioned is directly brought into 
question by the runs anomaly which was closely examined in 
reference 13. This property, defined in Chapter II. section 
A., informally states that under the symbolic dynamics, every 
possible sequence of zeros and ones is possible or realizable. 
If a pseudorandom number generator has a favorable runs 
property, all possible n-tuples (sequences of zeros and ones 
of length n) must not only be possible but their occurrences 
must be balanced. The unbalanced count of particular n-tuples 
1s precisely the runs anomaly pointed out by Forre as the 
fatal flaw of her pseudorandom number generator. 

The runs property of a binary sequence 1s tested by 
Seumeing the occurrences of the 2" different possible n- 
tuples. For example, in Heyman’s study, typical binary 
sequences of length 10% were used to test 2-tuples on a 
computer. There are 2‘=4 2-tuples: {0,0},{0,1},{1,0},{1,1}. 
The program ONTO.M (Appendix A) counts the four different 27 
tuples in a binary sequence. Figure 39 shows a bar graph 
Pee Gives the count of the four 2-tuples for a typical 


binary sequence of length 5 x 10°. 


55 


wv 
F) 
v 
a=) 
= 
5) 
a 


Z ee 3 
Shifted Decimal Equivalent 





Figure 39 Incidence of 2-tuples in typical Hénon 
generated binary sequence 


Notice that the horizontal axis is actually the decimal 
equivalent of the n-tuple shifted by one from zero. That is, 


instead of the following 2-tuple to decimal correspondence: 


2-tuple decimal 
00 ~ 0 where {0 0] * [2 1])* =e 
O01 - a where {0 1] *9(2°1)> sa 
10 - 2 where [1 0] * [2 1]° == 
alas - 3 where [11] * [2 1]° =a 
we use 
00 = a 
on = 2 
10 = 3 
AE = 4. 


We will use this shifted correspondence for all n-tuples. The 
unbalanced or uneven bins of 2-tuples in Figure 39 
demonstrates the runs anomaly. The depth of this anomaly is 


shown in Figures 40-42 which correspond to 3,4 and 5-tuple 


51S) 


counts respectively for typical binary sequences of length 5 
par 0. Modifications of ONTO.M were used to produce the 


associated data for finding the count of 2-tuples to i7- 


tuples. 
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Figure 40 Incidence of 3-tuples in typical Henon 
generated binary sequence 
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Figure 41 Incidence of 4-tuples in typical Hénon 
generated binary sequence 
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Figure 42 Incidence of 5-tuples in typical Hénon 
generated binary sequence 


Notice that in Figure 41 the bin corresponding to the decimal 
13 1s empty. The bin corresponding to the decimal 7 is also 
empty but in tests using binary sequences of length 10° Ehis 
onligueals) ielenz . The 13 bin, however, remains empty in the 10° 
length test. 

It was believed [Ref.13] that as sequences of greater and 
greater length were tested, the missing sequences would be 
found although the runs property, namely the balance of the 
bins, probably would not improve. Under memory constraints of 
the runs property test uSing a Sparc station 2, the lemGmamen 
the longest possible testable binary sequence was roughly 10°. 
Using initial conditions from all four quadrants in separate 
tests with binary sequences of length 10° no 4-tuple was ever 
found to correspond to “the decimal The decimal 13 
corresponds to the binary 4=tuple {1,1,7070}). Tters thereto 


expected that, when the 5-tuple runs are tested, those 5- 


So: 


@memes Which contain the 4-tuple {1,1,0,0} (see Table 2) 


should be empty. 


TABLE 2 5-TUPLES CONTAINING {1,1,0,0} AND DECIMAL EQUIVALENTS 








Pee LesSswithea{l»l,O» 04 decimal equivalent 


ie WO 07 0)) 2s | 
5017.0 seis) 26 
Ope, Tee ee 13 


{ dels, Oey 





This 1S indeed the case; however, these are not the only 


subsequences of length 5 unrealized in a sequence of length 
O°. Mable 3 lists those additional 5-tuples that are not 


realized and their corresponding decimal values. 
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TABLE 3 5-TUPLES AND THEIR DECIMAL EQUIVALENT 


= = - 


5-tuple decimal equivalent 


{09 07 OFC O 7 
(OF 07 Ue) 


(iy OG 


ss A tt = - oe — 


Cire 28 


In fact, the property was tested for n-tuples from n = 2 





to 17 (see Appendix C for 6 through 16-tuple runs) and an 
increasing percentage of unrealized sequences occurs. Table 
4 shows the number of unrealized sequences and the total 


number of possible n-tuples for a particular n. 
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TABLE 4 UNREALIZED SEQUENCES AND TOTAL POSSIBLE n-TUPLES 
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Figure 43 shows how the percentage of unrealizable sequences 


increases as the n-tuple length increases. 
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Figure 43 Percent of empty n-tuple bins vs. n 


B. ANALYTIC PROOF OF UNREALIZABLE FOUR-TUPLE 

Based on the previous data alone, the presumed 
homeomorphism cannot be disqualified from having the onto 
property. An analytically substantiated example of a wholly 
unrealizable sequence, however, would suffice. Let us 
investigate the first suspected unrealizable sequence, 
eee ee Os From the runs tests we see that the sequence 
{l,1,0,1} is possible although {1,1,0,90} 1s not. Fig@eeue 
shows the portion out of 5000 points on the attractor which 


gives the sequence {1,1,0}. 
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Exclusion Area of Points which Never Give Sequence 1100 





Figure 44 Points from a sequence of 5000 that give the 
sequence {1,1,0} 


These are also the same points which give the sequence 
{1,1,0,1}, since none gives the sequence {1,1,0,0}. hess 
apparent that these points are localized in a particular area. 


Figures 45-47 show the subsequent iterates of these points. 


Se 





Figure 45 First iterate of points that give the sequence 
3 aD) 





Figure 46 Second iterate of points that give the 
sequence {1,1,0} 
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Figure 47 Third iterate of points that give the sequence 
me, 0} 


It 1S clear that all these points under iteration correspond 


to the sequence {1,1,0,1}. 
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subquadrilateral or subquad). 
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Subquad, Quadrilateral and Split of. 4098 
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Figure 48 Subquadrilateral placed around those points 
that give the sequence {1,1,0} 


By iterating the points belonging to this subquad using the 
program FIND110.M in Appendix A (see Figures 49-51) we see 
computer-generated evidence that no points within the subquad 


correspond to the binary sequence {1,1,0,0}. 
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Iterated Subquad.Quadnilaseral and Split of .1098 





Figure 49 First iterate of points comprising 
subquadrilateral 


Iterated Subquad.Quadniateral and Split of .4098 





Figure 50 Second iterate of points comprising 
subquadrilateral 


67 


Iterated Sabquad.Quadrilateral and Split of 4098 





Figure 51 Third iterate of points comprising 
subquadrilateral 


This exclusion zone does plainly suggest that the sequence 
{1,1,0,0} 1s unrealizable. Using a binary Hénon sequence of 
length 10°, approximately 5115 points lie in the subaquacwern 
exclusion zone. This suggests that an orbit has a probability 
density of approximately 5.115% in the subquad. Since the 
Hénon map has been shown numerically to have topological 
transitivity, we expect that under reverse iteration of each 
point in the subquad there will be a preimage which will lie 
in the subquad. Therefore, this means that all points pass 
through the subquad and that no points which pass through this 


window give (2) 076). 
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Let us look more closely now at the points that would give 
meee, 0} imstead of those that would not. If we can show 
that the set of points that can possibly give {1,1,0,0} is the 
empty set, then we will have proved that this proposed 
symbolic dynamics does not possess the onto property. Recall 


the Henon map introduced in Chapter I: 


1 - 1.4x,° + y, 


a 
I 


ve = en. 


Meeeeoilows that: 

3G engl oe ee > coe Cp ce (4.1) 
Strerirst goal is to find the solution set that corresponds to 
ree}. In order for a sequence of (x,y) values to 
Sereees pond tO {1,1,9} the x values {xX,.;, X,, Xn.,} MUSt obey the 


following inequalities: 


Doane oor i (4.2) 

“), na AMed (4 ° 3 

See Se (a4) 
gene x,..=.4098, the dynamic median rounded to four 
Peemericant digits (previously discussed). Inequality (4.2) 


can be substituted in equation (4.1) as follows: 


ees aa) oo eee og (Aas) 


n+] 


oe) 


where the constant 1.12294 is eogneea to five signifies. 
digits as an underestimate in order to make our proof valid. 
Moreover, in all such cases we round in the appropriate a 
conservative) direction. In Figure 52 we see via (4.5)@@mac 


the solution set corresponds to the region above the inverted 


parabola, in the x,x,.,-plane. 





Figure 52 Solution set of points that could give the 


sequence {1,1,0}=umder <2], restrict aen 
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Figure 53 shows the possible solution set under restrictions 


fee (4.3) and (4.4) 


| 
| 
| 
| 





Figure 53 Solution set of points that give the sequence 
meee) under x,, restrictions 


By inequality (4.4), x,., has a maximum value less than or 
PemaumecO X,.,. Lf we let x,., equal exactly x,,, 1n inequality 
(4.5) then x, 1s further restricted to a value greater than 
Bea 1. Thus, the solution set corresponding to equation 


(4.1) 1s limited to the region shown in Figure 54 and reviewed 


as follows: 


X,-1 : region above the inverted parabola : 1 
Ce LOM SUCK that «.-> ./1371 eu: 
Pe ee GeglOn Such enae x; = .4096 220 


fas 





Figure 54 Solution set of points that JEner ace the 
sequence {1,1,0} 


By equation (4.0)" 26 follows scehiaeee 


PS 
r 


ne Lat Sr ee 
then 
5 Sea eee ae (4.6) 
To achieve the sequence {1,1,0,0}, x,,, must be less than or 
equal to .4098. However, 1t 1s clear that the largest x,,, 
value in equation (4.6) occurs at a minimum value for x, and 
a maximum value for x,,.,, which lies in the solution reqiommag 


Figure 54. The minimum x, value and maximum x value are 


nel 


respectively .71371 and .4098. Using these optimum values in 


equation (4.6), 


V2 


eee ee fA Seely) = be a 098)“ 
=e 57 9003. 
The symbol (+) .1ndicates that we round to the number of 


Significant digits shown and that we round in the appropriate 
(conservative) direction as previously defined. ies) be 
Merve (not surprisingly) corresponds to the point with the 
Meamum x value in Figure 51, that is, the third iterate of 
the subquad region which only gives the sequence {1,1,0,1}. 
Of course, this minimum x,., value exceeds .4098, and 
therefore, the sequence {1,1,0,0} 1s not realizable. Thus, we 
have found a counterexample to the proposed homeomorphism. 
Furthermore, 1t 1S shown in Appendix D that the sequence 


moe, 0,0} 1s also analytically unrealizable. 
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V. CONCLUSIONS 


A. RUNS ANOMALY-FATAL FACTOR OR LIMITING FACTOR? 

We have seen that the classical Hénon map exhibits the 
attributes normally accepted as characterizing Ween. 
dynamical systems. These properties include sensitive 
dependence on initial conditions, topological transitivity and 
a dense set of periodic points. Despite this, however, we 
have provided ample numerical evidence and rigorous analytical 
proof that the proposed*symbolic dynam@es scheme h: TR Seen 
generating pseudorandom binary sequences 1s: 

1) not a homeomorphism for the proposed symbolic dynamics 
scheme since it 1s not onto. (Because the Hénon attractor 1s 
a subset of TR we can also conclude that h:H, ~ 2 is noe 
ele eer 

2) highly restricted in its viability as a pseudGrangem 
number generator. 

We have shown that not all binary sequences are generated with 
equal frequency. In particular we have shown that certain 
sequences are not realizable and that others are very sparsely 
attained. These facts support the observations and 
conclusions of Réjane Forré which suggest that this scheme is 
unsuitable as a reliable means of generating pseudorandom 


numbers. 
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The evidence suggests that the runs property severely 
limits the potential of the scheme for widespread practical 
use. It 1S shown in reference 18 that subsequent elements of 
the binary sequence may actually be predicted with 70-80% 
em@euracy by an artificial neural network (compared to 503% 
Seemracy for a coin flip), most likely due to the severity of 
@@emeuns anomaly. However, because of its Simplicity and 
nonlinearity there may exist some applications for which the 
scheme would be well-suited. It 1s our belief that the 
general idea of using a chaotic discrete dynamical system to 
generate pseudorandom binary sequences, Nowever, has merit and 


deserves additional study. 


B. FUTURE WORK FOR AN IMPROVED SYMBOLIC DYNAMICS 

Our numerical results reflect the basic structure of the 
attractor. Despite the lack of a homeomorphism the structural 
mature of the attractor is apparent in the binary output of 
the symbolic dynamics scheme. AlEHOUGH Che “attractor Ais 
accepted as possessing chaotic attributes, we believe that due 
Memea= Structure of this “chaotic driver" (the classical Hénon 
map) pseudorandomness is not fully realized in this scheme. 
We believe that a sequence of more than four zeros 1s not 
possible because the left fixed point is both repelling anc 
Me@teeom the attractor. Contrarily, it has been shown that a 
sequence of up to 23 ones 1s possible [Ref.16]. We conjecture 


(as suggested by computational evidence) that the sequences of 
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ones are possible because the right quadrilateral fixed point 
is a saddle pemne- Specifically, the sequence of ones is 
possible due to the attracting axis of the saddle point. This 
suggests that the symbolic dynamics scheme could be effective 
if there were such a saddle point on the attractor on both 
sides of the dynamic median Xycn. 

Since the location and nature of the fixed points@ie mean. 
Henon map depend on the (a,b) parameters of the map, 
[Ref.1:p.170] there may exist parameter sets which provide 
this structure. The chaotic bands present in the bifurcation 
diagrams suggest there exist many other (a,b) pairs that may 
give even more complicated dynamics. Proving the existence of 
(a,b) parameters that give us these desired characteristics 1s 
only the firseceeo- It is also required that a bounded 
attractor exist for the associated parameters. Furthermore, 
Forre had one eae of being provided a numerically 
calculated trapping region from which to take the pseudoranmaem 
number generator key. This region, if 1t exists at al ewenme 
have to be recalculated for the new (a,b) pair. 

If it could be proved that there does not exist an (a,b) 
pair which gives rise to saddle points that possess the 
previously mentioned attributes, then it may be possipiemme 
use another map which corresponds to a different dat tiagman 
with these qualities. However, the major asset of the Hénon 
recurrence iS its utter simplicity which translates tauiame 


fast generation of pseudorandom sequences. 
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Basing the symbolic dynamics scheme on the x value of each 
meetore 15 Simple but it is not clear that it is the most 
effective. If we are willing to accept a more complicated 
S@ememe we Could base our sp/iit on a different linear median or 
a nonlinear median. The goal in choosing a given median would 
MeweorOring parity to the runs property for the system of dual 
saddle points while retaining the other properties previously 
mentioned. A more complicated scheme might generate binary 
sequences more slowly, but the loss of speed may be warranted 
to improve the runs property. 

It may be argued that despite an improved system there 
will still exist binary sequences that are unrealizable. 
However, the severity of the associated runs anomaly could be 
diminished should the shortest of these unrealizable sequences 
be of sufficient length. That 1s, we have seen that an 
unrealizable 4-tuple was catastrophic in this case because the 
problem of unrealizability translated to every subsequent n- 
tuple. In fact, we must anticipate that as described in 
Chapter IV even more n-tuples will be unrealizable than 
expected. The degradation of the system thus 1s directly 
related to the length (L) of the first unrealizable 
sequence(s). Since a pseudorandom number generator 1s judged 
by more than just the runs property criterion (see conclusion 
Ref.13) a system which possesses a "large" L-value may still 


provide an acceptable pseudorandom number generator. 


(af. 


APPENDIX A: MAIN PROGRAMS 
NOTE: GRDCOMP PROGRAMS REQUIRE THE HENREAL PROGRAM TO RUN 


function [xmat, ymat] = grdcompl(sp,W,Z) 

This function takes a point(W,Z)from the RT quad,initiates 
Giac! cadc 

% field based on a certain spacing, then iterates those grid 
field points 

* that match the binary string of (W,Z) [use({(.2),.1) feGrine we 
As long as 

$ the strings match, points are iterated using the Henon 
recurrence.) OMe ena 

%$ (those) points that completely match are finally plotted 
according, towrne 


of oP 


o) 


NO 


% following scheme: 

% lem Griginal Grider relem pene. 

% tt Or LiGisde ae 2 

% go 1terated field points with matching binary string 
abenin iteration 

S wt iterated (W,Z) at nth iteration 

% Ops original xvec,yvec point(s) which matches binary 
$ string of W,Z. We also show how W,Z and 


xvec (indices), yvec(indices) “walk" to the Henon attractor at 
each iteration (with the attractor on screen). 

o=linspaiece (—-1es, 1 esc, 500 

Ss==. 106220. wee O- 

u=lingepacei@s, 3 27a245 Siar 

Voor 64 * mee 4. Oy ho- 

G=linspace(—1205, be245 S00. 

Ql Se) 2) tof Says al 

e=linspace(-1.06,-1.33,500); 

fom. 407 *c -4, 11 oe 


Solit = eeaesise 
a = [-. 1082 =1; 3. 6451] 256752). 
b = [-.2760;-4.6718;-.3344]; 
GC = (Spine: sprees cn 
ee |— 6: spree |= 
[x,.y¥].-= ~aenweal (70s. — iG es 
ximene = [Split splice). yline = [=1982 1.3278 
lenc = lengrmre. lend = seme bl c). 
xVeC = = zeros(1,lenc* Jemma 


zeros(1,lenc*lend) ; yvec 
LE mint(a*|we2Z)) =o == Ue | oes eae 
disp initial valwmewms nietean the rent quad: | 


geleubaghs 
end 
eos le 
FOR 4p) Bee Lene 
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meget = sl: lend 


Kea Kee 1s 
mame (a (CD) 7a ij < 5b) == 0 
Bepee (kK) = Welle avec ik) = d(x); 
else 
magee(K) =e Use yectk), = —-10: 
end 
end 
end 
fee = reshape (xvec, lenc*lend, 1); 
yvec = reshape(yvec,lenc*lend,1); 
m = find(xvec ~= -10 | yvec ~= -10); 
pee = xXVEC(m); yvec = yvec(m); 
oe = Py aigeler. Ve Vee 
WwW = WwW; 
22 are 
peewee |- 1.32 1.32 -.6 .5]); 
meeaees = 1:length(xx); ee 2 ee 
1 =1; 
meec = lLength(xx); 
ee 2haLces: 
leng = length(xx) % 
iemetalization of gq 
Mmmm oovec, vec, b.’,W,2, 'r*’,xline,yline, ’w-’);hold on; 
Peete, Ss, W.' U,V, 'W.',;6@,f,'’wW.'’,09g,n,'’w.’); Rold on: 


meee Grid field and W,Z’);pause;clg;hold off; 

ame Jength(q) > 1; 

eeero([-1:32 1:32 -.6 .5]); 
) 


x0 = yy + ones(size(yy)) - 1.4*xx.%2; 
pte. SO * XX; 
ie = 272 + | - 1.4*ww.%*2Z; 
mee = .3* ww: 
Peewo > Solit 
Caen senior oplit | x0 == cp lase) 
else 
eet = L£inci<0 = splat): 
end; 
Deelengch qteme) == 0, break; 
else q = qtemp; 
end 
mem it: 1: $ counts iterates where at 
least one point matches 
leng = length(q) 
qvec = [qvec,length(q)];% plot of how # with same binary 
sequence 
ma@ereess = 1ndices(q) ; %$ decreases each time through loop 
pee = x0 (Cd); $ preserves to next iterate those 
matching values of 
We = YO a); $ x0 and yQ 
ww = w0; 
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27 ee 


pause; 
plot (xx,yy,’g.',w0,z0, wt’ ,xline,yline man |) holawear 
plot(o,S,'’W.' U,V, We. 4 ey bees 9 Cpe ee se olee ane: 

title(’Movement of points in grid with same binary 

sequence’ ) 
pause> cla; nolagerr. 

end 

ax1s; 

vec = (1/2). ((1 lenge aiGgves 


piot((1:length(qvec))-1,qvec, ‘r’, (1: length(aqvec) )-1/qveqgiiaa 
Veer. ).] 
title(’Decrease in # of grid points with same binary 
sequence vs 1/2%n’) 
pause;clg; 
$ In this way, we need the index key ‘indices’ to tell to 
Wilehe Grid pemme weme 
$ Surviving iterate corresponds with respect to the original 
xvec and yvec. 

axis ( (=) 3 2a ome 


plot (xline,yline, ‘w',x,y,’b.’,0,S, ‘w.’,u,V, 'W. 72,1, ee 
i Woe Oe sone, 
plot (xvec(indices) ,yvec(indices),‘'gx’,W,2Z, 'r*’) ; ReOweaae 
title(’W,4 and Ene grid spoeinme warm Ssaomemoumam, 
sequence’) ;pause; 
he Fe or re 
Knew = xvec (indices) 
ynew = yvec(indices); % now we show a plot of how W,Z and 


Ene peint 
% that generates the same binary sequence 
xnew, ynew 
(@) 


%$ walk around the attractor 
axis. ( [ea ea Sia 


plot (xline,yline, ‘w',x,y,'b.’,0,S,'wW.’,U,V, We. we, Le 
pO oe 
homer, 
XnNewtl = XNewW; 
ynewt = ynew; 
WE = WwW? 
ial Cae — a 
FOr) n= 167i 
KOC, S50 el er ewer ee en ewig 
Albom Ps sn daleyiic : 
WOE c eO ele ee 2 t eeee 
Z0C = eo ie 
Wea WUE; 
ZC ea eee 
XnNewe  “=scUiar 
ynewt = yOt; 


plot (MOt,VOt/ "Gx ,wOt) 20 ena 
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title(’Movement of pt chosen on HA and matching sequence pt 
until diverge’); 

meme, plot (x0t,y0t,’ix’,w0t,z0t,’'1i*'); 

end 

aera Ort; 
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function [xmat, ymat] = grdcomp2(sp,W,Z) 

This function takes a (W,Z) from the left quad,initiates a 
Wade Cis Geb em! 

based on a certain spacing, then iterates those grid field 
POLlnes mae 
% match the binary string of (W,Z) [use ({(-1,-.25) for merase 
As long as the 
%$ strings match, points are iterated using the Henon 
recurrence. On hema 
% (those) points that completely match are finally plotted 
according Ger ene 


oe oO 


oP ,Q 


% tollowing sehneme: 

% ioe Orl ginal “Guid! Ziela woernes 

% es Original (W,Z) 

% ofe) iterated field points with matching binary Seema 
ab Hehe eerat 10, 

% W+ iterated (W,Z) at nth iteration 

% opr original xvec,yvec point(s) which matches binary 
Stile sol eee 

% We then plot the iterates of W,Z and the point that most 


closely matches 

$ its binary sequence. 

O= linspace(—fy see 2 Seu 
= Oe Ot o. 

UW ieems pace (i 32 e 45a) Cee 

Via) Cs TIL 

G=linspace( —l0e7 e245 7500) 

Bewilose 40-2 oo 

e=linspace(—-1.06,-1. 225500). 
f=—3.407*e-4.- 1 lie 


Sle = eee 
a= [344094 1; = 106351 33567 ae 
Doe = 4 IS oO = ea 
GC. =p es 2 eo So late 
olga lee sie) pl at 
[syeey | = benreal (7 50, =1707..-2 >. 
enc = lhewctrn (cc). lend = bength va: 
wine. Seemed “Sjodeante |e viane = |-).32 1032) 
xvec = zeros(1,lenc*lend) ; yvec = zeros(1,lenc* lends 
VE mint(a* [wrz2l) S95) == 05 i ee rae 
disp(’initial value is not in left quadr’) 
Per Uiem. 
end 
k=O 
fOr p=) lene 
LOL? stele Ke el 
k=k+1; 
1f°max((a* (6(p)5 dita ne) a 
MVEC(K)S=NC( Oi we Se le eae ia 
else 
xvec(k) ="=107 yvec(h) =) 50, 


end 
end 


reshape (xvec, lenc*lend, 1); makes a column vector 
reshape (yvec,lenc*lend,1); 
Eame(xvec=-= =10 | yvec== -10) ; 
paec (In) ; VVeC== VVec (mM); 
coe XVeC; Vea eC: 
WW = W; PT hon Saad 
peeeewe 32 1.32 —.6..5)}); 
indices ne Telia 4 < ) an |e ae oe ee 
ae ; 
qvec IGhals pelau..<)as 
| inca Cces > % 
[metallization of g 
meng = length(xx) 
melot (xvec, usta; 2). ae Zz efor ee eee ho ldvion: 
Meter4O,S,'’wW.",Uu,V, ' WwW.’ Pal aed Gy eens) Ow Lr. 
eee be(’ Grid eles and W, Z’ eee Cig; holed orr- 
Peele lengthiq) > 1; 
epees( (—-l232 1. ov Ono) 
( 


< 

D 

QO 
nou Ww dl 


St ee 


mae = YY + ones(size(yy)) - 1.4*xx.°2; 
hae, = .3*XX;> 
feme= 22 + | - 1.4*ww. 2; 
wee .3*ww; 
Merwe =. Split 
‘Meeomow= site <0 > "splat |) x0 == spdat); 
else 
A@eCiipn = rome x0. < sp lie) 
end; 
Moslengriaiqremp) == 0, break; 
else q = qtemp; 
end 
m= 1 + L: 
moe = length(q) 
ec = (qvec, length(q) ]; 
moermees = i1ndices(q); 
eae =" <0 (q) > %$ preserves to next iterate those 
matching values of 
yy = yO(q); %$ x0 and yO 
wi = WO; 
wae = ZU; 
pause; 
Meb@rixx,yy, o.'’,wO0,20,’w+’,xline,yline, 'w-');hold on; 
Peueierre. We 4 tl, V7, We pe. t,’wWwe’,0,h,'’w.’); nNold off; 


mee |’ Movement of points in grid with same binary 
sequence’); 
pause;clg 
end 
ax1s; 
feereme= (1 /2).7>((1.:length(qvec))-1); 
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plot((1l:length(qvec))-1l,qvec, ‘yr’, (1: length (qvec) ) -l)qVeesanue 
VieC Veo owon 
title(’Decrease in # of grid pts with same binary sequence 
as leer IMGs) 
pause;clg 
% In this way, we need the index key ‘indices’ to tell to 
WALGh- Grid perms tae 
% surviving literate corresponds with respect to the original 
xvec and yvec. 

axis |(=bo32) 82 oes ee 
plot (0,8, 'W.*,U,V, W.', eC, ly We Po, We ene eo me 
plot (xvec(indices) ,yvec(indices) , 'gx’,W,Z,'‘'xr*’,x,y, 3 
e,yline,‘w’); 
title(‘’wW,Z and the grid point with the same binary 


sequence’) ;pause; 
hole rete: 
xnew = xvec(indices); 


ynew = yvec(indices) ; 
axis ( (=1.327. e323 =. 6 oer 
plot (xline,yline,’w’,x,y,'bD.',0,S,'W.',U,V, WwW.’ , 2, 2 


xnewt = xXnew; 
ynewt = ynew; 
Wir awe 
7 he Tage — 
Com 7 6) 
“OL = Lby0 = 1. 4*xnewte 2 erie 
VUER= 2.3 newe: 
WOG> Slee Oe = awl Jeter 
Z.\) Gs ee Gee 
WE = woe. 
Zi. =a 
xNewt = 50 


YVOewe = 520 
DIOt(XOUGFVUL, Ox’ ,WOU 20b) bee lcmear 
title(’Movement of W,Z and matching sequence pt until 
divergence’) ;pause; 
DLOL (| x0G5 VU ewe eee 
end 
lake woes ie 1s 
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% function [xmat, ymat) = grdcomp3(sp,W,Z) 

% This function takes a point(W,Z)from the RT quad and 
1terates the point n 

% times in order to ensure the point 1s on the Henon 
attractor(ensuring 

meoae the nth iterate is in the left quadrant). It 
initiates a grid field 

% based on a certain spacing, then iterates those grid field 
eetees Chat 

* match the binary string of (wnew,znew) [use (-1]1,-.25) for 
Vee As long as 

% the strings match, points are iterated using the Henon 
recurrence. Only that 

% (those) points that completely match are finally plotted 


e@@eercing to the 


%$ following scheme: 

% liom Seimei t1eld points 

% ial Srigqimal (Wiz) 

S go iterated field points with matching binary string 
at nth iteration 

% w+ iterated (W,Z) at nth iteration 

% gx Original xvec,yvec point(s) which matches binary 
Sjeemmg of WwW, Zz 

% We also show how W,Z and xvec(indices),yvec(indices) 


walk" to the Henon 
attractor at each iteration (with the attractor on 
screen). 

@oeeemsmoace(-].33,1.32,500); 
Pees *O + 2/76; 

Wewmememace(1.32,1.245,500); 

Peet -4~.6718; 

seumeacoace(-1.06,1.245,500); 

Peewee oO sy 3344; 

e=linspace(-1.06,-1.33,500); 

memos *e -4.1119; 

seme = .4098; 


ove 


peepee OSS -1;-3.64 1;-.1562 1}; 
meee = 2 /760;-4.6718;-.3344] ; 
@meemereplit:sp:1.32]; 
Smee O7Sp: .5); 
meme = henreal(750,-1.0,-.25); 
Peeeme = {split split]; Vine =F (=P 52 2]; 
eames = lengthi(c); lend = length(dqd); 
Pevee) = zeros(1,lenc*lend) ; vyvec = zeros(l,lenc*lend) ; 
keme=  O: 
Pome = |1:lenc 
meme re=. |: lend 
k=k4+l1 ; 
mamiiceoe a | C(pD) side) )] i =< b) == 0 
Mae Gu) s=) Cio) -iyvec (kk). =(dir); 
else 
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xvec(k) = -10> =-yvec@e) = ies 


end 
end 
xvec 
yvec 
m 
xVeEC 
yvec 
Wace 
Zae eis as 
FOr aa ea we Theis 
ensures wattr,zattr 1s 
wattr(j+1) = 1.0 - 1.4*wattr(j)*2 + zattr())> 4emiee 
abe ct Ou 
Zaeen a+) lee eo 1 Wa tea a 
1f j3 > 20 & wattr(j+1) > split, break; 
end 


reshape (xvec, lenc*lend,1) ; 

reshape (yvec,lenc*lend,1); 

find(xvec~= -10 | yvec~= -10); 

xvec (m) ; 

yvec (im) 

Wi % This preserves the values of W,Z 


end 

wnew = wattr(j+1); $ preserves value on attractor 
WARE, Zatti 

ZW — ww Zo wees alee a 
Walter = whew; 

zattr = 
xXx = xVec- $ preserves the values of xvec & yvec 


Vee, 2.32) So Gee 

i lenge i (se)5: + [1 2. 3 Saeed 
i 
ence (xx): 

G] inGgieces: 
leng length (xx) % initialization Oma 

plot (xvec,yvec,’b.’ ,wattr,zattr, ‘r*’,xline,yline, ‘w= eee 
Om, 

plot(0,S,’wW.' ,Uy Ve We’, Cr EyY We’ 7G, ae We Onno 
title(’Grid field and point chosen on attractor’ ) ; pauses 
while length(q) > 1; 


avec 


axie{ (—1326) 32 sa ome se 
ORS" Yyy + Ones | Siaetiaw as = a ee 
Wa = a ee 
wO = zgnew + 1 - 1.4*wnew.“2Z; 
ZZ). = ee ie 
1 WO eee 
qtemp = find(x0 = splice woolen 
else 
GiemMmp \=ehilne 0 2 sp lena 
end; 
1f length(qtemp) == 0,break; 
else q = qtemp; 
end 
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P= ot |S %$ counts iterates where at least 
one point matches 
meee = (Qvec,length(q)];% plot of how # with same binary 
sequence 


leng = length(q) 
me@mecs = 1ndices(q); $ decreases each time through loop 
Pe =) x0 (Gq); % preserves to next iterate those 
matching values of 
bare = YO (q); % x0 and yO 
wnew = wQ; 
@aow = zQ; 
pause; 
Meme, YY, ‘oa.’ ,w0,Zz0, ‘'w+',xline,yvline, ’w-’);hold on; 


Mme >, Wo’ 7 U,V,’ wW.',e,f,'’w.',d,n,'’w.'); 
title(’Moveméent of points in grid with same binary 
sequence’) 

m@emise; Ccla;hold off; 


end 
feeeen | l 2 3 4)) Tal Ss 
meer (1/2). re enlencente. cen = lee: 
mbeee lt: length(qvec))-l1,qvec,’r’, (l:length(qvec) )-1,qvec(1) * 
meen b’); 


Eitle( ‘Decrease in # of grid points with same binary 
sequence vs 1/2%*n’) 
DWause;clg; 
$ In this way, we need the index key ‘indices’ to tell to 
Pome orid point the 
paeeeyiving iterate corresponds with respect to the original 
xvec and yvec. 

peo |-1.32 1.32 -.6 .5]); 
feimeoryne,yline,’w',x,y,'b.’,0,sS,'w.',u,v,'’w.’,e,f,'w.’,g,n 
meee) Old on; 
meeemecyec (indices) ,yvec(indices),’gx’,wattr,zattr, ’r*’);hold 
om; 
Mellen Wattr,zattr & grid pt with same binary 
sequence’) ;pause; 


nomea, Off; 
pilewe= <vec (indices) ; 
ynew = yvec(indices); % now we show a plot of how 


wattr,zattr and the point 
%$ that generates the same binary sequence 
xnew, ynew 
% walk around the attractor 
eee - 1.32 1.32 -.6 .5)); 
Memecline vyline,'’w',x,yY,'bD.’,0,S,'w.',u,v,'’w.',e,f,'w.',g,n 
ew ) 


mela On; 

xnewt = xnew; 
ynewt = ynew; 
Poeeret = WaLlttr; 
memere = Zactr; 
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for n= eee 

xOt = 1.0°- 1. 4*xnewre 2. neve. 
Vow = 3° xnewe, 
wO0t = 1,0 =9)-4*Wwaterwun 2) Zac wie 
20%. (=o Noe oe 

Wattle eae wo ee 

Zatert, ={22 0 ee 

XNewt = xX0ET 

ynewt = yOt; 


plot (x0tyy Ob Ox 2s vii 0 epee a 
title(’Movement of pt chosen on HA and matching sequence pt 
until diverge’); 
pause: “plot (x0 thy 0 1 eee eee 
end 
jaliaplicl Ges 2 
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% function [xmat, ymat] = grdcomp4(sp,W, Z) 

% ONLY DIFF BTWN 3&4 IS THAT 4 WEEDS QUADR,3 USES WHOLE 
RECTANGLE. 

peegas function takes a point(W,Z)from the quadr and 
iterates the point n 

%$ times in order to ensure the point is on the Henon 
attractor(ensuring 

meeemat the nth iterate 1s in the left quadrant). It 
Memerdates a grid field 

% based on a certain spacing, then iterates those grid field 
femmes that 

$ match the binary string of (wnewW,Znew) [use (-1,-.25) for 
VWeeoee AS long as 

% the strings match, points are iterated using the Henon 
recurrence. Only that 

% (those) points that completely match are finally plotted 


Bee@eraing to the 


% following scheme: 

% D. Srrgqitial gqurdelreld points 

% ft Serger (WZ) 

% go iterated field points with matching binary string 
Semiten 1teration 

ra wt iterated (W,Z) at nth iteration 

% op. Original xvec,yvec point(s) which matches binary 


Seemng of W,Z 
Eemimat long 
% We also show how W,Z and xvec(indices),yvec(indices) 
"walk" to the Henon 
mete eractor at each iteration (with the attractor on 
screen). 7 
@mispace(-1.33,1.32,500); 
Seas *O + .276; 
Meeeerspace(] .32,1.245,500); 
meee. -4.6718; 
@eteimepace(-1.06,1.245,500) ; 
ete o~*g -.3344; 
eeeamsopace(-1.06,-1.33,500); 
Beno. 40/*e - 4.1119; 


seer = .4098; 
ewes 2074 1; -.1083 -1; -.1562 1); 
foe 4 1119; -.2760; 2 ols 
pee). 32:Sp:split]; 
emem—.6:Ssp:.5]; 
mmr mete [w;z)) > b) == 0 | w > split 
Saistot dmleral Value 2s not in lett quadr”) 
return; 
end 
ee y) = Nnenreal(750,-1.0,-.25); 
mimic = (Split split]; yee see | 
enc = length(c) ; kenidy= lengtintd) 
pavec = zeros(l, lenc* lend); yvec = zeros(1,lenc*lend) ; 
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% This next section puts everything outside of the quad to 
(= Os z— nes) 


if “Max ( (a. “Set ea as een oe — 
xVeEC(K) ="e€(p) =yveo Ge) =e aa )s 
else 
xvec (kK age) -10; Staged (ko. = a 
end 
end 
end 
XVeC = reshape (xvec, tenes lends 
yvec = reshape(yvec, enc lena aa 
m = find(xvec == —-DOMP (veer. - je 
xVeCC = <Vee (mie, 
VVeCC =) y Vee (mi 
WAEED = - % This preserves the values of W,Z 
Zee =e 
EOne) =e % 
ensures wattr,zattr 1s 
wattr(j+1) = 1.0 - 1.4*wattr(j)*2 + zattr(j); 3 Onwieee 
attiraceou 
Zauer (a+!) = Vo wattria.. 
1£ 7 > 20. 8 wWaterig+!) = spolie moreak 
end 
end 
wnew = wattr(j+1) ; 
ZROW = Seale isa) le: 
wattr = wnew; 
ZAGEr = Z2new- 
So ‘=o cee %$ preserves the values of 
xvec & yYveCc 
YY — sa Vee, 
ane [=1.32 1375-6 
Indices. Je lengen is. 2 [253 44 
ae 
Gqvec = lengiat-<aoi- 
q = indices; 
leng = length(xx) 


plot (xvec,yvec,’b.’,wattr,zattr, 'r*’,xline,vline, Ww jeeeeee 
on 

plot(o,s,'wW.',U,v) WwW.’ 767 £, Ww. 7G, ny ee oltclern ar 
title(’Grid field and point chosen on attractor’ });pause7eumer 
while length(q) > 1; 


axrSigp- 132 eZ Oreo 

x0 = Yy +eOnes (size) eel ae 
YU = <2 to. 

wO = znew + 1 - 1.4*wnew.%2; 

ZzQ0 = .3*wmew; 
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PET WwO Ss spileee 


ae Toa ct Sole Wie) == split): 
else 
Qeemo Sneenc |e i Sela & |; 
end; 
Pee length (qeemp) == 0, break: 
else q = qtemp; 
end; 
= 1+ 1 % counts iterates where at least 


one polnt matches 
qvec = [qvec,length(q)];% plot of how # with same binary 
sequence 


Meng = length(q) 
ieeme@ees = indices(q)  ; % decreases each time through loop 
fece= <0 (cq): $ preserves to next iterate those 
matching values of 
a= VYU(a) -; oS. <0 eane. 0 
wnew = w0Q; 
anew = zQ0; 
pause; 


meow xs yy, '9.',w0,z20, w+’ ,xline,yline, ‘w’);hold on; 
manele Ss, WwW. fu, WwW. ey lp we’, dg, n, ow." )+ 
title(’Movement of points in grid with same binary 


sequence’); 

meuce> Cig; hold off; 
end 
ax1s 


epee 3 4)); axis; 
ee 1/2).*((l:length(qvec))-1); 
Meerlengren(qvec))-l,qvec,’ xr’, (1l:length(qvec) ) -1,qvec({1)* 


title(’Decrease in # of grid points with same binary 
Sequence vs 1/2“n’) 

pase; cla; 

$ In this way, we need the index key ‘indices’ to tell to 
Weel Grid point the 

$ surviving iterate corresponds with respect to the original 
payee) and yvec. 


peepee 1.32 1.32 -.6 .5)); 


meee ane, yline, 'w’,x,y,'b.’,0,8,’w.’,u,Vv,’'w.’,e,f,’w.’,g,n 
mee. )s; hold on; 
mmeienxvec (indices) ,yvec (indices), '’gx’,wattr,zattr,’r*’); 
ager on; 
Meee Wwattr,zattr & grid pt with same binary 
sequence’) ;pause; 
Mmewa oft; 
Pewee = xvec(indices) ; 
ynew = yvec(indices); % now we show a plot of how 


hemeer,zZattr and the point 


2) 


$ that generates the same binary sequence 
xnew, ynew 


° 


Sewoaleearouna the.attractor 


gi 


e 


= 
oie 
Ges 
ame) 
Ks 
Siok Bl 
Hl 


Zee ee 
POlg th 2= al 
<Ot— 
VOr 
WOE 
Zc 
Wale ei 
2 came 1 te 
xnNewCt = 
ynewt = 
plow 


OZ, Lao Ober ean ie es 


hold em 
xnew; 


aL 
Vy, Db. jxllane, vlane, we fo, see ve eee 
; 


= ynew; 


Wael. 

Zot cies 

1 - l 
1 WO vss Ls 4 Wir cee ey eo Wie- 
no  SHewe: 
20) S=i4 twat eee eee nee 
yo Welle oe 

Vides 

Zita 

Setar 

eeiey 

ROG, VOC, OX" Gwe 2 Oe) 


ne. 


',e,f, wees 


title(‘Movement of pt chosen on HA and matching sequem@ema. 
Giga sions ence fae) 


pause;plo 
end 
NnolesenEr- 


C(x0t VOR. ix swe eee 
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function [xmat, ymat] = proofl(spcrs,spfn,incr) 

meeeo FUNCTION’ GOES BACK TO **GRDCOMPI1** AND USES IT OVER 
iMeniy SP VALUES. 

This function takes a point(W,Z)from the RT quad,initiates 
sa quad grid field based on a certain spacing, then iterates 
S$those grid field points that match the binary string of 
ae) (use{.1,.1) for now. As long as the strings match, 
points are iterated using the Henon recurrence. Only that 
(those)points that completely match are finally considered. 


OO oO 


WP oO 


of oc oX 


a=. 2; 
oe .1; 
ooieee = .4098; 
mepeeo = SDfn:incr:spcrs, %$ MAJOR OUTER LOOP 
foe] SD 
pee i ees 3-15-3264 1;-.1562 1); 
mel. 2700: =4.6/18 >=. 3344); 
Beso Ltt Sp ls, 32]; 
fo ls 62 Sp? 5]; 
Memes= length(c); lend = length(d); 
pevee@e— zeros(1,lenc*lend) ; Vvee = Zemoc( lo benc™ lend); 
mame (a*(Wwez)])- =S= b) == 0 | w < split 
Sisco initial value is mot in the raght quadr”) 
feeeurn > 
end 
income et a Oe 
meee 1:lenc 
meneame = |: lend 
k=k+1; 
mua <a Peto). sa1r) |) -< by ==) 0 
MCenin) SC (py ay vee (kK) = air); 
else 
Meyeo (ky oO A Vvec (kis s= = 10; 
end 
end 
end 
mvec = reshape(xvec, lenc*lend,1); 
yvec = reshape(yvec,lenc*lend,1) ; 
mma f£ind(xvec ~= -10 | yvec ~= -10); 
eee  xVeC(m); yvec = yvec{m); 
oe = XVeC : VY = yvec; 
WW =e ZZ =a 
Mmeemmees = |: length{xx) ; Se oe og: 
32" as Jails 
ereec = length (xx) ; 
=e 1ndi ces: 
feng = length(xx); 


while length(q) > 1; 


mue= Vy + Ones (size(yy))-1.4*xx.72; 
Pace 2 3 x: 

eee zeae = 1 Aww. Oe: 

BG. = a8 WW: 


32 


Lf wO0 ss sepia 


qtemp = find(x0) > spire | Ue —eoe ae 
else 
Gtemp = fimGtecs isola yy 
end; 
1@f Jbengen( qecme  S==" 0" brea. 
else q = qtemp; 
end 
ee cles Le $ counts iterates where at 


least one point matches 
Leng =. Lenore nia) . 


qvec = [qvec, length(q)];% plot of how # with same bamene 
sequence 
IMG Gesa—wl We he@ecenre %$ decreases each time through loop 
yO ep <0 eC) $ preserves to next iterate those 
matching values of 
Wa 20 Ce $ x0 and yO 
ww = wQ; 
ZZ Tera 
end 


° 


% In this way, we need the index key ‘indices’ to tell to 

WAiCh Gr idepol tmeieme 

$ surviving iterate corresponds with respect to the original 

xvec and yvec. 

xnew = xvec(indices) ; 

ynew = yvec(indices); % now we show a plot of how W,Z and 
S$the point that generates the same 

binary sequence xnew,ynew walk around the attractor 


xnewt = xnew; 
ynewt = ynew; 
WE = we 
Ze Sys 
EO pile: ele eee 
KOC =) ta0 = 1. 44snewre. 2 ee ewe, 
VOC =e. slewior 
WOE = 120s 1 .4*We72 Zee 
20 @ =) 23 ee 
We =-Ww0e: 
Zt KS Ue 
Snewe snes 
Views | Jee: 
end 
nvee += |mvecem|- 
spvec = [spvec,sp]; 
clear xvec;clear yvec; 
end % END OF MAJOR LOOP 
axils; 


pLOE(Ssbv eG aVee ean 
title(*’SP VS. ITERATIONS UNTIL BPilVSrGenee a. 
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function [xmat, ymat] = proof2(spcrs,spfn,incr) 
Mero FUNCTION GOES BACK TO **GRDCOMP2** AND USES IT OVER 
Mey SP VALUES. 

femme LUunction takes a (W,Z) from the left quad,initiates a 
Mea grid field 

based on a certain spacing, then iterates those grid field 
Seats that 

Meee the binary string of (W,Z) [use (-1,=-.25) for now]. 
As long as the 

%$ strings match, points are iterated using the Henon 
recurrence. Only that 

% (those) points that completely match are finally 
considered. 


oe oP 


0° TH of OQ oe 


eoe= —-1.0; 
Z, =o ee 
peeme = .4098; 
moepemeeo = SDfn:incr:spcrs, %$ MAJOR OUTER LOOP 
eo) = Sp 
ember 407 127-—.1083 —1l--~-,1562 LL); 
fomem(—4.1119;- =.2760; -.3344]; 
moet - 1 3275p -splxt |; 
al -.o:Ssp:.5]; 
Meme = length(c); lend = length(dqd) ; 
xvec = zeros(l1,lend*lenc); VWvecC.= zeros (17 lene vend); 
meampmiaiia*(wez)y > by == 0 | w > split 
Gees initial value as mot in the Vert quadr’) 
return; 
end 
a 
aero = |:lenc 
moe t= 1:lend 
=e: 
immencocta” (C(p)-dir))}) = b) == 0 
Pee Coe GD ay vec (Kk) <= Ging); 
else 
mvec(K) =~-10; yvec(k), = :-10; 
end 
end 
end 
xvec = reshape(xvec,lenc*lend,1); $makes a column vector 
yvec = reshape(yvec,lenc*lend,1); 
Mme Lind(xvec~= -10 | yvec~= -10); 
Seer = xvec(m); yvec = yvec(m); 
noe =Exvec; VY = Vaec: 
WW Sey ae =e, 
meemecs = 1:length(xx) ; oe ee eee. Was 
1 = 1; 
@vec = length (xx) ; 
Qe = andaces; $ initialization of q 
mene = length (xx) ; 


mele lengthigq) > 1; 


25 


x0 = yy + ones(size(yy ))))— 234 oe 
yO = «303284, 
wOo = Zz + 1 eae 
2 =) ee 
Li wO 23s lie 
qtemp =] find (0 => spike ee S== ce eee 
else 
Geemp = “kine c0s = semen. 
end; 
1f length(qtemp) <== 0,break; 
else q = qtemp; 
end 
je a eae 
Peng = hemo ieke i. 
Givec 9=0 (Gee wen gtn ie an 
Lid Cece =a riGumees cy). 
pe aa a Ooi) %$ preserves to next iterate those 
matching values of 
ave All sce 6 x0 and y0 
ww = w0; 
eg ar as 8 


end 
% In this way, we need the index key ‘indices’ to tell to 
Wilh Ged deren tate me 


$ Surviving iterate corresponds with respect to the original 
xvec and yvec. 


xnew = xvec(indices); 
ynew = yvec(indices) ; 
xnewt = xnew; 
ynewt = ynew; 
WE v=" Ww. 
ZG = Ze 
fOr n=, laa =a | 
KOC = 1. 0.<) 124*xmewt.. -29e newie, 
Vue. = YWasxnewe: 
wOt = 1,0 = 1.4*wer?2? ec 
ZUG. = eee ee 
yi = Chie 
Ze, = Santee. 
Mew  =se0it 
ynewt = y0Ot; 
end 
nvec =“fmwec. ml 
SPVEC =] Fseveo ook 
clear xvec;clear yvec; 
end % END TO MAJOR OUTER LOOP 


DLOEISPVeEC aieG. 2 yaaa 
title(’SP VS. ITERATIONS UNTIL DIVERGENCE’); 
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$ function [xmat, ymat] = proof3(spcrs,spfn, incr) 

pape FAN USES (W,Z)=(-1,-.25). ONE ATTR PT IS COMPARED 
OVER MANY SP VALUES 

Peeeiiey DIFF BIWN 3&4 IS THAT 4 WEEDS QUADR,3 USES WHOLE 
RECTANGLE. 

@etnis function takes a point(W,Z) from the quadr and 
iterates the point n 

%$ times in order to ensure the point is on the Henon 
attractor(ensuring 

% that the nth iterate is in the left quadrant). It 


mnitiates a grid field 

based on a certain spacing, then iterates those grid field 
mBeimes that 

%$ match the binary string of (wnew,znew) [use (-1,-.25) for 
W,2]. As long: as 

* the strings match, points are iterated using the Henon 
Me@mrrence. Only that 

° 


“eeeose) points that completely match are finally 
considered. 


- 


ow 


oe . 9; 
wee . 1; 
pees = . 4098; 
Weletr = WwW; $ This preserves the values of W,Z 
Soeer = Z; 
memes; = 1:40 % 
P@euees Wwattr,zattr is 
Pomme |) 6 v= 1.0 = 1 a*wattr(7)°2 + zattr(s); % on the 
attractor 
maeer ijt!) = .3*wattr(z}; 
marge 20 & Wattrigtli > split). break. 
end 
end 
wnew = wattr(j+l); 
eieweme = Zattr( +1); 
eepemom = Spin:incr:spcrs, % MAJOR OUTER LOOP 
ee = sp 
wattr = wnew; 
zattr = znew; 
emee|-.1083 -1; -3.64 1; -.1562 1]; 
meee — .2760:-4.6718; -.3344] : 
me sp lat:sp:1.32); 
mee |-.60:S5p:.5); 
ogee =—= lengthi(c)  ; meme = cLengin (di. 
Pegeees— zeros(l,lenc*lend) ; VWVec = 72er0s (),bene “lend).; 
Mamma iw;zi) > b) == 0 | w <= split 
Meso Initial value 1s notin right quadr’) 
Teale Opel 6 Or 
end 
oO 
Meme = 1:lenc 
mem yr = 1 lend 


oY) 


k=k+1; 
$ This next section puts everything outside of the quad 
( : i. 


to (-10,-10) 
1f mMesaiia * ,[c(b) dir) |) -aeie— 
XVEC (Kh) aC yi CeiwG) = vener) | 
else 
xvec(k) = -10; yvec(k) = -10; 
end 
end 
end 
XVEC = reshape (xVvec bene; lends 
YVEC = reshape (yVvec, Heme Vena. 
m= find (xvec.—= —L0 VS veer 7 Ole 
xVEC = XVEC (Mj) VeCc = YVee my. 
Oe = eee YY Pe a 1S 
indices = deVengen (is) S ible 24 ee 
se — erly: 
GVeCe =—= lenctim 2.0 
q = indices; 
Jeng = "ebeng El. e: 
while length(q) > 1; 
x0 = yy + ones (s1Z2 (a7 eis * xx. 2; 
V0) =. ae 
w0Q = znew + 1 - 1.4*wnew.%2; 
2.0 =e AS wie. 
16 WO) S= "Splat 
gqtemp =9f1nd (x0 =) sp late vine) ee ee 
else 
GEemo = £1nd (<0 Sp lenieae 
end; 
1f length(qtemp) == 0,break; 
else gq = qtemp; 
end; 
a ele $ counts iterates where at least 
one point matches 
gvee = 'laqvec, Length. (qi) ie % plot of how # with same 
binary sequence 
Teng = length (ai 
InGitess > = saancmces (qs % decreases each time through loop 
co. = ee): %$ preserves to next iterate those 
matching values of 
Vag Oo s $ x0 and y0 
wnew = wQ; 
ZnNeW = zZ0 ; 
end 


° 


6 In this way, we need the index key ‘indices’ to tell to 
which grid point the 

6 SUYViving iterate corresponds with respect to the original 
xvec and yvec. 


Hew = eVec i nemecs 


ge 


ynew = yvec(indices); %* now we show a plot of how 
Peter, Zattr and the point 

6 that generates the same binary sequence 
xnew, ynew 


oe 


walk around the attractor 


xnewt = xnew; 
ynewt = ynew; 
oeeet = Waltr; 
Meer t = zattr; 
age = 13 - 1 
tc 1.0 - 1.4*xnewt.*2 + ynewt; 
moe = . xnewt ; 
Mien ele — 1.4*wattrte.°2 + zattret; 
eee = .5*wattrt; 
Wattrt = w(t; 
meee = ZOt; 
Pg@ewo = xXO0t; 
Powe = YOt; 
end 
meee = (nvec,n); 
fevec = {[spvec,sp]; 
@teagmexvec; clear yvec; 
end > END TO MAGOR LOGE 
aAX1S; 


mega (Spvec,nvec, ‘'w*'); 
mete oF VS. ITERATIONS UNTIL DIVERGENCE’); 


* function [xmat, ymat] = proof4(spcrs,spfn, incr) 

* THIS FXN USES (W,Z)=(-1,-.25). ONE ATlTRee) 1S COMeawa 
OVER MANY SP VALUES 

% ONLY DIFF BTWN 3&4 IS THAT 4 WEEDS QUADR,3 USES WHOls 
RECTANGEE 

This function takes a point iW,Z) from thesquacdwecme 
1terates the point n 

$ times in order to ensure the point is on the Henon 
attractor (ensuraund 

% that the nth iterate is in the left quadrant). It 
Initiates a grid field 

% based on a certain spacing, then iterates those grid field 
points that 

$ match the binary string of (wnew,znew) [use (-1, -. 25 0a 
W eclity AS. LONG vas 

%* the strings match, points are iterated using the Henem 
Yeecurrence. Cn enae 


* (those) points that completely matcenm are tina, 
Cons Vdered. 


W = —l1SG. 
ZS =e 
S OG ae ee 
Wak eae %$ This preserves the values of W,Z 
Za Wie — ze, 
EOr eae eo % ensures wattr,zattr 1s 
wattr(j+1) = 1.0 - 1.4*wattr(j)7*2 + zattri(j); 2 omen 
ack racuen 
Zactrigel) “= yee water (4) 
if 3 > 20 & wattr(j+1) < split, break; 
end 
end 
wnew = wattr(j+1); 
NeW = GzZaetr (ae le. 
[Or Sp = seins ner seems %$ MAJOR OUTER LOOP 
SOs = SD 
wattr = wnew; 
ZaAcGwe= Zhen 
ans. | seedy .d. 1 rere se =) ae eae 
Be= [eae a2 C0 eo ale 
ope re cle? se) ismen| - 
elee=. “are See oe 
lene =: TenGgente). lend = length(d) ; 
XVECGU= Zenos (l,Lene Lema yvec = zeros(1,lenc*lengs 
LE mini ta* (wsz)) S ey ==) 0 iu cele 
disp(‘initial value 1S net aWnetes mouse. | 
return; 
end 
ke Oe 
[Or Dp =e. lene 
[Omen melee) Clad 
9 eno ae 


g89.0 


eethis next section puts everything outside of the quad 
EemeeanO , -10) 
1f max((a * a. ) ] 
Reet) <= .e (p)s ec 
else 
Ree (6) -=— LO eyvec(k) =.-10; 
end 


se Oy) -=- 0 
(OS = ist aen 


end 

end 
ot @ 
yvec 
m 
-eVeC 

ond 

indices 
cl 
qvec 


reshape (xvec,lenc*lend,1); 
reshape (yvec,lenc*lend,1) ; 
ema (xyec ~= -10 | yvec ~= 
xviecim);  yvec yvec(m); 
xVeEC; var yvec; 
ies bengeh (xx) > lw. eek, 4s Gee 
ie, 
ema t ti (sx ).- 
g IMaLees ; 
leng STIG) 
mwomee length(q) > 1; 
ba et Ones isizetyy)) = 1:4*kx.°2;> 
mo x; 
maew + 1 = 1 .4*wnew. 2; 
.3*wnew; 
Peawos= Solit . 
Cmeme S7 enews = so lit. "|KO 22] split)s 
else 
@eemp = find (x0 ysplit). 
end; 
Pembengeitqtemp) == 0, break: 
elece CG = °9qQtemp; 
end; 
eo 1 + 1: %$ counts iterates where at least 
one point matches 
@eec = [aqvec, length({q)]; % plot of how # with same 
binary sequence 
Heng = 
indices = indices % decreases each time through loop 
ream= XO (cq); % preserves to next iterate those 
matching values of 
VY MOC) ; % x0 and yO 
wnew w0; 
znew cad Oe 
end 
6 In this way, we need the index key ‘indices’ to tell to 
which grid point the 
 SUrViving iterate corresponds with respect to the original 
Pee and yvec. 
xnew = xvec(indices); 


oer atte oral 
HA 


< 
oO 
How ow oat 


uO 1 


ynew = yvVeciemd1 Gees now we show a plot of howwattr, 
zattr and the point that generates 
the same binary sequence xnew,ynew 


walk around the attractor 


of oP oO AO 


xnewt = xnew; 
ynewt = ynew; 
WAEULE = “Waeer, 
Vis SUG) @) Cue — MEAs EC 
FOr sie — ole: ee 
xOt-= 10 = Pea *2newe. 2a weve 
Vt = 2233 Slevin, 
wOt =s1-70 = 1. 4*wattrie Zao ese 
ZU oes Swale ie 
Wattrt = w0e-; 
Zaria 2G. 
Knew = 205. 
VMewe “i =) vy vies 
end 
NveG we lh Vee oa. 
Spvec- = Ulsovec spi 
Clear xVeG. cleanly vee 
end % END TO MAJOR LOOP 
anus. 


DEOE (SUveC avec we 
titie(‘’SP VS. ITERATIONS UNTIL Di VERGENGES © 


G2 


% function [xmat, ymat] = prfl(spcrs,spfn, incr) 
peas TunNction takes thirty-one points (W,Z) from the RT 
quad,initiates a 


co) 


peace grid field based on a certain Spacing, then iterates 
those grid field 


° 


mme@esncs that match the binary string of (W,Z) {use(.1,.1) 
Ber now. AS long as ; 
% the strings match, points are iterated using the Henon 
recurrence. Only that 

(2) 


% (those) points that completely match are finally 
considered. 


aoe = O- 
eam = .2:.02:1.1 $ MAJOR OUTER LOOP 
ee ee 
mm 11 + '1 
peepee = .4098; 
meee = Spin:incr:spcrs, % MINOR OUTER LOOP 
ee = Sp 
papeee > 083 -lies. 64 1;=-.1562 1]; 
meee. 27600;-4.6718;—.3344]); 
mee Soe Sp. 1 32): 
eos | -.6°Sp:. 5]; 
Meme |= lLengthni(c); lene v= kenge (d).; 
Peer zeros(!,lenc*lend) ; yvec = zeros(l,lenc*lend) ; 
eminem (a- (wez)]) > Db) == 0° |) ow <= Split 
SmcO atitial value 12S not in the right quadr’ ) 
ie ebb eoak 
end 
aro: 
mere = |:lenc 
mer ry = |1:lend 
k=k+1; 
MINCE atch LC) Ga tineacs) == 0 
Muee we)! = ielp)s so voet kK) [=< (r} ; 
else 
reel | =O Veeck) 25 = 20- 
end 
end 
end 
fage@e= reshape(xvec,lenc*lend,1); 
yvec = reshape(yvec,lenc*lend,1); 
farina (xvec ~= -10 | yvec ~= -10); 
eee XVCC(M); Yyvec = yvec(m); 
YOK =e VveCc: VY = yvec; 
WW = Ane Z2 — ee 
maces = 1:length(xx); Be A ee eae ee al 
ao = i) 
eevee = lLength(xx) ; 
Ge = mo vees : 
meag = length({(s~x) ; 


m0 3 


while length (qm 1: 
x0 = yy * Ones (stzely7) ae eae 


57.) = ena ere: 
w0 = zz + 1 - 1.4*wwl 
20) = ae 
Blin ee eh ouleage 
qtemp. = fand(x0 = splie i) x0) ==. een 
else 
Sleeiiic amas iglollio<di ez > Siellkavie i" 
end; 
1£ lenGgthiigqtemp) == sR bweac: 
else q = qtemp; 
end 
ee %$ counts iterates where at 


least one point matches 
leng = length(q) ; 


qvec = [qvec,length(q)];%3 plot of how # with same binary 
sequence 
LGAUCeS == anctee sta) %$ decreases each time through loop 
oer > JUN onie: $ preserves to next iterate those 
matching values of 
Vy = setae = <0 andy? 
ww = wQ; 
Zot a en: 
end 


° 


%$ In this way, we need the index key ‘indices’ to tell to 
which grid point the 
$ surviving iterate corresponds with respect to the original 
xvec and yvec. 
xnew = xvec(indices) ; 
ynew = yvec(indices); % now we show a plot of how W,Z and 
the point 

$ that generates the same binary sequence 
xnew, ynew 


9° 


% walk around the attractor 


xnewt = xXnew; 
ynewt = ynew; 
We = WwW 
ZZ. =a 
fie ee ole eee | 
KOC = 1.0 - 1.4*snewe. 20S news 
VUE = 23° <Mewe: 
WOEs= 12.0 = lea) eee eee 
ZG eee Wise 
WANN adore tees 
Zi — ee 
gla iewy — >is 
VOCwW Ces —— ee 
end 
nvecl © = Jnyceitaa.- 
Sspovecl = i[spvect ce 


104 


clear xvec;clear yvec; 


end % END MINOR LOOP 
meec2 = |[nvec2,nvecl]:; 

spvec2 = [{spvec2,spvecl]; 

clear xvec; clear yvec; 

end % END MAJOR LOOP 


axis([0 max(spvec2) 0 max(nvec2)]j); 
perou(spvec2 ,nvec2, ‘w.’); 
mere’ SP VS. ITERATIONS UNTIL DIVERGENCE’); 


igs 


function [xmat, ymat] = 


of oO 


MANY SP VALUES. 


3 (AND OVER THIRTY -ONGn. 2 cLlUES, 
%$ This function takes a (W,Z) 

quadHe rudd. = tere 

%$ based on a certain spacing, then 
points that 

% match the binary String of ee 


As long as the 


° 


recurrence won, sevar 


% (those) points that completely match are finally 
considered. 
sheers (C) 
[Or gw = =250 24 %$ MAJOR OUTER LOOP 
m= = Oe 
ica leee eA 
30) Jot ee COS Fe) 
FOr. Spy] (Spin Ines suemce % MINOR OUTER LOOP 
Sp = Sp 
a= (3.407215 =.1083—) aoa ils 
Db = (4.1 a 2 io eo oa 
oP —NIIe Ae sicit Sia Jaie 
C=. Oe See aie 
lene = lLengmne : lend =" Vengen (ad); 
“MVEC = zerosiil, lend. lence) |: yvec = zeros(1, lenc* lenge 
LE Man (tat wee 1 aa Oy ee so lee 
disp(’initial value is not in the left quadr’) 
return; : 
end 
K=O. 
EC) aes Che — lea als 
fOE % =] 1) lend 
esis alii 
iPemax ((a* (Cp) dias) alo ee 
sevec(k) =e Vp) avec ie) = eae: 
else 
RVEC (kK) t= "Ue avec i. =~ LU 
end 
end 
end 
xvec = reshape(xvec,lenc*lend,1); $makes a column vector 
yvec = reshape(yvec,lenc*lend,1) ; 
mM =TbaGiGevec-= -l0 0 7ee-— ter 
XVEC = Vee ny; yVeC s= 7 Cen e, 
oe =xVec- VY = Vavee 
WW = We Ze = <7 
indaces = leona ela. 2 lee ee 
iy = ils 
qQvec = lenge aec. 
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prf2(spcers,spfn,incr) 
THIS FUNCTION GOES BACK TO **GRDGOMP2<~- 


AND USES IT Gia 


from the left quad,initiates a 


iterates those grid field 


[use (-1,-.25) for mea 


% strings match, points are iterated using the Henon 


oe 


S= 2@ndices: 
Mereialization of q 
Meng = length(xx) ; 
Woeee lengthi(q) > 1; 


me] yy + ones(size(yy)) - 1.4*xx.%2; 
Pathe 45 "XX; 
wee 22 + 1 = 1 .4*ww. 2; 
moe 3 3 * WW; 
enw > Splat 
ecto m— meri Kure Solit. | «Oh == split )-- 
else 
Geeme v= vr Pieee.8 e— “Sore 
end; 
1f£ length(qtemp) == 0,break; 
else q = qtemp; 
end 
=e 1 + 1: 
Wemg = length(q); 
@yvec = [qvec,length(q) ]; 
maemces = 1ndices(q) ; 
poe <0 (cc); $ preserves to next iterate those 
matching values of , 
yy = yO(q); $ x0 and yO 
ww = wO0; 
mame 7,(): 


end 


° 


$ In this way, we need the index key ‘indices’ to tell to 
Piemeneoriad point the 
° 


$ surviving iterate corresponds with respect to the original 
vec and yvec. 


p@liew = xvec(indices); 
ynew = yvec(indices) ; 
xnewt = xnew; 
ynewt = ynew; 
we. = WwW; 
Zee = 2: 
mete >i - 1 
miime= 1.0 - 1.4*xnewt.°2 + ynewt; 
fee = .35*xnewt; 
ieee 1,0 — 1.4*wt.°2 + zt; 
mee = ,3%*wt; 
me = wOt; 
mie = 70 t 
fmewte = xOt; 
ynewt = yOt; 
end 
mec = (nvecl,n}; 
Spvecl = {spvecl,sp]; 
clear xvec;clear yvec; 
end > END “POeMINOR OUTER LOOP 


OF, 


nvec2 [nvec2,nvecl]; 
spvec2 {spvec2,spvecl1]; | 
end % END TO MAJOR OUTER LOOP 
ax1s([0 max (Sspvec2) (0 mmaseimveczia | ; 
plot (Spvec2 nvec2 eae 
title(’SP VS> 1TERATIONS UNTIG PiLVERGENCE 12 


seeks) 


% function [xmat, ymat] = prf3(spcrs,spfn,incr) 

fee > FAN USES THIRTY-ONE W,Z VALUES. ONE ATTR PT IS 
COMPARED OVER MANY SP VALUES 

SeONLY DIFF BTWN 3&4 IS THAT 4 USES LEFT QUADR , 3 USES 
RIGHT QUADR. 

fm@ihes function takes a pointS (W,Z) from the quadr and 
lterates the point n 

$ times in order to ensure the point is on the Henon 
attractor (ensuring 

$ that the nth iterate 1s in the left quadrant). It 
Maeetaces a grid field 

% based on a certain spacing, then iterates those grid field 
points that 

Pe@atecha the binary string of (wnew,znew) [use (-1,-.25) for 
W,Z}. As long as 

$ the strings match, points are iterated using the Henon 
recurrence. Only that 

% (those) points that completely match are finally 


@emcidered. 


m= 0; 
Mme = .5:.02:1.1 $ MAJOR OUTER LOOP 
Zz ae 
[mie = 11 + 1 
Seeeec = .4098; 
wast = WwW; $ This preserves the values of W,Z 
wee = Z; 
more) = 1:40 % ensures wattr,zattr 1s 
weueerig+l) = 1.0 - 1.4*wattr(j)°2 + zattr(j); % on the 
emaaclOr 
meee 7t+l) = .3*wattr(7) ; 
Memeo 20 & Watt GE (ee) JS: Soli tye break: 
end 
end 
wnew = wattr(j3+1); 
Bigeweee= Zattr(j+1); 
Memes = Spin:incr:spcrs, % MINOR OUTER LOOP 
See= sp 
wattr = wnew; 
Zee = ZNew; 
meme es =-1; -3.64 1+ —.1562 1); 
f@e= |—.2760;-4.6718;-.3344]; 
Men cer tt -sp:l.32); 
ema |. 6 7Sp:.5]; 
femenr= !ength(c); tence] -engen.(d) 
Pome e= zeros(l], lenc*lend) ; Vijieece= Zeros. lence lend), 
memmun((a*(w:z}) > b) == 0 | w < split 
Gesell Inteial valwe tis mot in right quadr’} 
: return; 
end 
i 


mee fm = 1:lenc 


Lew 


for r = bavene 
k=k ee . 
$ This next section puts everything outside of the quad 
Cali 10) 
1£ max((a =) [Clone d(T) | eee — 0 
xvectk) =9eN—p) eck = aia 
else 
XVeEC(K)) = -=10- 3 Vy 7ee ti eee ee 
end 


CO 


end 

end 
xVeC 
yvec 
m 
xVeEC 

< 

indices 
a 
qvec 


reshape (xvec,lenc*lend,1); 
reshape (yvec,lenc*lend,1); 
find(xvec ~= -10 | yvec ~= -10); 
xvec(m); yvec yvec(m) ; 
xVee: VY YVeC; 
ie emiG ean ex > |1 2 eee eee) 
It, 
enolofe lal (see "- 
q indices; 
leng IBeugtefe leo 
while length(q) > 1; 
x0 Vy + ones(sizeliyy)) = 2le4 7s. 2 
yO ee 
w0 znew + 1 - 1.4*wnew.%2; 
ZU .3*wnew; 
if WO m= =o So leints 
qremp = find (30 =—splwe (2:0 —_-) ool. 
else 
QCeMpe= Lind <0r< Spies 
end; 
if length(qtemp) “== 0; break 
else G-="Gtemnp-. 
end; 
ge ee Ale $ counts iterates where at least 
one point matches 
Gqvec = ihevec, fengt nh (ee $ plot of how # with same 
binary sequence 
beng = Jength( a) 
indices = indices(q); % decreases each time through loop 
POE = Ocean. % preserves to next iterate those 
matching values of 
YY VOs Co $ x0 and yO 
wnew wO0; 
znew ZO 
end 
% In this way, we need the index key ‘indices’ to tell to 
Maulelal cresel jelecinic ela 
$ surviving iterate corresponds with respect to the oriageaaes 
xvec and yvec. 
xnew = xvec(indices) ; 


| 


ee 


ynew = yvec(indices); % now we show a plot of how 
Petr,zattr and the point 

% that generates the same binary sequence 
xnew, ynew 


ov? 


walk around the attractor 


xnewt = xnew; 
ynewt = ynew; 
Pee etet = WaCtr; 
meepetk = Zacttr; 
aemeem = 1:1 - 1 
Pie = 1.0 - 1 .4*xnewt.*2 + ynewt; 
eae = .3*xXnewt; 
Poe = 220 - 1.4*wattrt.°“2 + zattre; 
wet = , 3S" wattrec; 
wate rt = wOt; 
@eeteiae = z2O0t; 
Pmewe = xQ0t; 
pmewt = yOt; 
end % END TO MINOR LOOP 
meee! = ([nveci,n]); 
peewee! = [spvecl,sp]; 
emeeeie XVEC; Clear yvec; 
end % END TO MAJOR LOOP 
meee = ([nvecz2,nvecl]; 
spvec2 = [spvec2,spvecl]; 
end 


axis([0 max(spvec2) 0 max(nvec2)]); 
Meemispvecz,nvecz,’.’); 
meee.’ SP VS. ITERATIONS UNTIL DIVERGENCE’ ) ; 


icra 


% function [xmat, ymat] = prf4(spcrs,spfn, incr) 

* THIS FAN USES THIRTY-ONE W,2Z VALUES Oh Ea eee 
COMPARED OVER MANY SP VALUES 

% ONLY DIFF BTWN 3&4 IS THAT 4 USES THE LEFT QUADR, Ja ue 
THE RIGHT OUARRe 

% This function takes thirty points(W,Z)from the quad and 
1terates the point n 

%$ times in order to ensure the point is on the Henon 
attractor(ensuring 

% that the nth iterate 1s in the letesquadrane ere 
initiates a grid field 

%$ based on a certain spacing, then iterates those gri@m@igeud 
DOINES Cnet 

% match the binary string of (wnew,znew) [use (-1, —-.2Z5)uee 
WZ). Skene as 

% the strings match, points are iterated using the Henon 
PeCurrence.. Gly sume 

% (those) points that completely match are finally 
considered. 


aces ae 
LOr Wt] >a.o0 es % MAJOR OUTER LOOP 
a eel 
BRO Meslay Secll 
Splits] 2420ce. 
Wat En = <we 6 This preserves the values of W,Z 
Gee et == 72s 
one wg = acl s %$ ensures 
Swattr,zattr is 
wattr({j+1l) = 1.0 - 1.4*wattr(])*2 + zattr()); seg 
4att raGeew 
ZACEr (eee eee oe ee ee 
2£ 3 > e2URe wattrijel) = solu» meee 
end 
end 
wnew = wattr(j+l1); 
ZEW Se eae er lanl 
EOr Sp = Splemeiney-soens, %$ MINOR OUTER LOOP 
So = se 
wattr = wnew; 
2attm = Z2hew 
ae= (324074 1; ~-.0esee esol eae 
Bees [=4e 11S: = 276? taneeA A. 
ei (ledes2 2S soles 
Git") [Saree spe" Se 
Leng = seme lh Geil. lend =—-lengtihnicie 
xvec = zeros(1,lénc*lena) ; yvec = zeros(1,lenc*lend® 
1f mint{la* [wz S95) == 70 cole 
Gisp(’initial value is “mot 19 lereeeiad. |) 
LeeWUen: 
end 


pie 


femeo = 1:lenc 
moe Fr = j:lend 
k=k+l1> 
$ This next section puts everything outside of the quad 
to (-10,-10) 
mee (a et oe te ecco) =" 0 
ee Wieomed al Vvectk) = d(r); 
else 
Ruece Ky = = 0" vveci ky) = =10- 
end 
end 
end 
xvec = reshape(xvec,lenc*lend,1); 
Sayec = reshape (yvec,lenc*lend,1); 
meer inda(xvyec ~= -10 | yvec ~= -10); 
xVeC = xvec(m); yvec = yvec({m) ; 
ond = XVEC; VY = yvec; 
maeees = 1:length (xx); eS ae ae 
1 = 1; 
Meee = length(xx) ;: 
@e- 2 indices; 
meme = length (xx) ; 
Peeme length(q) > 1; 
pa vy + Ones {sizel(yy)) - 1.4*xx.%2; 
wae = =. Ot KX? 
Tie= znew + 1 — 1.4*wnew.*2; 
2) = .3*wnew; 
Paw = Split 
Gqeewoe= a icioc’ = Splaiteie <0 SSeep lit )+ 
else 
Seco. — Piece or < Splat) 
end; 
tea lengua, auemp) == O0,break: 
else gq = qtemp; 
end; 
me 1 +4) % counts iterates where at least 


one point matches 
qvec [qvec, length(q) ]; 
binary sequence 


— 
= 


leng = length(q); 
icGgkieces = indices(q); % 
Pam = x01 CG) ; % 
matching values of 
yy = yO(q); % 
wnew = wQ; 
emew = z0-; 


end 
% In this way, we need the 
Meee Grid point the 


2) 


6 


xvec and yvec. 


Surviving iterate corresponds with 


° 
© 


plot of how # with same 


each time through loop 
to next iterate those 


decreases 
preserves 


xO and yO 
to tell to 


ice Key “indices * 


respect to the original 


Ls 


xnew = xvec(indices) ; 

ynew = yvec(indices); % now we show a plot of how 
waltr,Zzattr and the pormme | 

% that generates the same binary sequence 
xnew, ynew 


oe 


walk areunad the art rect om 


xnewt =. new 
ynewt = ynew; 
Watt bey = Wak eae 
Za Glee Sh: Zeilee tar 
for sn = 1 see 
xOt = 1.0 = 1. 4*xnewe Ze ever. 
Vito =" 23 ewe: 
wOt- = 1.0 =] 1.44 watt 3 Zech egies 
GWE. St Laas Wreieate Tae 
WalGne =<wlt- 
Zeabiibea 2 Ou 
OO =) ee — dice: 
VTLew teas ee 
end 
nvecl = (mv7eerrn | 
spvecl (= ifspveci spi. 
clear xvec; clear yvec; 
end % END TO MINOR DGGE 
Mveecz2 = [mvec2;nvecl|, 
Sevec2a = |[spvecz, sevechi: 


end % END TO MAJOR LOOP 

axis([0 max(spvec2) 0 max(nvec2)]); 
DIGE(SDVECZUnVees se We Je 

title(’SP VS. ITERATIONS UNTIE ot eeeo len ae 
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ONTO .M 

THIS PROGRAM REQUIRES THE PROGRAM HENON TO RUN PROPERLY 
[mars Program 1S a derivative of a program used originally 
in reference 13 called "runs.m". 

mors Drogram 1S run using the word-"onto” 

Been use "x = henon(#,0,0,)"; then use "onto"; 

for the appropriate n-tuple sought. 

memen (500000, 0,0); 


“ by SX o© ad a A of of Of 


= 4 
ee 
Peeevec2 = zeros(4,1); 
$decvec3 = Bares (a. a)"; 
@eevec4 = zeros(16,1); 
p@leevecS = zeros(32,1); 
saecvec6 = zeros(64,1); 
saecvec/ = zeros(128,1); 
$decvec8& = zeros(256,1); 
meleevec9 = zeros(512,1); 
$decvecl0O = zeros(1024,1); 
meeevecill = zeros(2048,1); 
$decvecl2 = zeros(4096,1); 
p@eecveci3 = zeros(8192,1); 
$decvecl4 = zeros(16384,1); 
$decvecl5 = zeros(32768,1); 
reece 16 = zeros (65536,1); 
meeeveci/ = zeros(131072,1); 
i= (> 
memebe min(decvec4) == 0 
ieee 
weeemmal = (2 1] * x(1:14+1); 
$decvec2 (decimal+1) = decvec2 (decimal+1) + 1; 


e@eocimal=(4 2 1)*x(1:1+2) ; 
$decvec3 (decimal+1) = decvec3(decimal+1) + 1; 


@eemmal=(8 4 2 1)]*x(1:1+3); 
@eevec4 (decimal+1) = decvec4(decimal+l) + 1; 


epoaeeimal=(16 8 4 2 1)]*x(1:1+4); 
peeeGyvecs(decimal+1) = decvecS(decimal+1) + 1; 


ee@eetmal=(32 16 8 4 2 1] *x(1:1+5); 
sdecvec6(decimal+1) = decvec6(decimal+1) + 1; 


peecimal=(64 32 16 8 4 2 1)*x(1:1+6); 
malecvec?(decimal+l) = decvec7 (decimal+l1) + 1; 


meeecimal=([128 64 32 16 8 4 2 1)*x(1:1+7) ; 
$decvec8 (decimal+1) = decvec&(decimal+1) + 1; 


1A 


$decimal=(256 128.7564 32706062 ee ec 
*decvec9 (decimal+1) = decvec9(decimals1) eer 


$decimal=[(512 256 128964 32°16 6342 ee 
$decvecl10(decimal+1) = decvecl0(decimal+1) + 1; 


decimal=[(1024 512 256 128 64 32 16 8 42 1/*x@ ee 
$decvecll(decimal+1) = decvecll1(decimal+1) + 1; 


$decimal=[{2048 1024 512 256 128 64 32 16 8 4 2 1)*x( 72a 
$decvecl2 (decimal+1) = decvecl12 (decimal+1) + 1; 


$decimal=[4096 2048 1024 S12 2565128964037 5o oe ee 
i pe a ae 
$decvecl3 (decimal+1) = decvecl13 (decimal+1) + 1; 


$v=(8192 4096 2048 1024 512 256 128 64 325632 oe 
$decimal=v*x(1:i+13); 
$decvecl4(decimal+1) = decvecl14(decimal+1) + 1; 


Sv=(16384 8192 4096 2048 1024 512 256 128 64 32 16 832 
2CeCima l= ae eee el) 
$decvec15(decimal+1) = decvecl15(decimal+1l) + 1; 


$v=(32768 16384 8192 4096 2048 1024 512 256 128 64 325 eee 
Cele lis: 

S0CGi mal =] 7 ae 

$decvecl16(decimal+1) = decvecl6(decimal+1) + 1; 


=(65536 32768 16384 8192 4096 2048 1024 512 256 lZegeae 
$decimal=[v 32 16 8 4 2 1] * x(1:1+16); 
$decvecl7 (decimal+1) = decvecl17(decimal+1) + 1; 


end 


$ 1 tells the sequence length at which all the n-tuples were 
Found: 


n %mn-tuple looking for 


1 % length of decvec where found 
length(x) % total length of original vector from Henge 


mr 


epumetion xx = find110(n,x0,v0) 

wes Program requires that plot(x110,y110,'.’) be used 
after the run. This program takes n points on the attractor 
and finds those that give the sequence 110. It then 


Siterates the set 3 times. Merely set the parameters below 
eee type "“findll10". 
m= 5000; 

mo = 0; 

yO = 0; 

Srmputs: 

% n = length of desired sequence 
% Peter a ba lx 

% Poren Initial y 

Bemeputs: 

% een lobia nary veelor 

mela On; 

oat) = x0; 
tae) = yO; 


pewmmsoace(-1.33,1.32,500); 
eee s*O + .276; 
meemasoace(].32,1.245,500); 
foe *1 - 4.6718; 
Sees pace(-1.06,1.245,500); 
Meeeesso~g - .3344; 
e—emesoace(-1.06,-1.33,500); 
maemo? *e - 4.1119; 


eee = .4098; 
peemee- 1-32 1.32 -.6 .5)]); 
feemem@ee— (Split split}; yline = [-1.32 1.32}; 


oer x line,yline, ’w-‘');hold on; 

meee, Ss, 'W.',uU,V,'wWw.',e,f,'’w.’,9g,n,'w.'); 
title(’Exclusion Area of Points which Never Give Sequence 
ieee O" } 


pause; 
ree) = zeros(n-1,1); % vectors are preallocated 
Pee, = zeros(n-1l,1); 
$recursSlve generation of points 
meme = 1:n: 
paweenreaHey (a) =F = [eAése( 1) 2: 
Vio = a * (1) 


° 


Seeonvert to binary 
mich) a= sso lit 

xe t= Oe 
else 


ae 


oe ele 
end 
end 


% Soe 
Geevecs 
1 


= ZOLOs Lo eae 

= 0; 

while min(decvec3) ==0 
i =i+41; 

decimal = 


L£ decimal==6- 


end 
end 


z Oke) ae © ciRIE Cea Ibs) Ae 


ite 


[a2 Ll ee eee 


pe IIS, 
eee 


[xl10,x( ge 
[yl10 yy Goes 


SUBQUAD 
This program shows how the four sides of the 
ubquadrilateral 


Py DL . 


oP oo M oP WK) oP of 


iemruin type "“subquad" 


fae inspace(-1.33,1.32,500) ; 
pee —.1083*xx + .2760; 

peer inspace(!.32,1.245,500) ; 
ies. 64*aa - 4.6718; 

weet anspace(-1.06,1.245,500) ; 
emer 1 533*cc -,. 3344; 
memeeeanspace (-1.06,-1.33,500); 
mms 5. 40/7*ee - 4.1119; 

Pee tanspace(.4099, .5714,500) ; 
poe. 1a*xX + ©2785; 

mee weeinspace(.5714, .5273,500); 
wee /O9"*a - . 8123; 

meen space (.5273, .4102,500) ; 
ome. 260/3*C + .26145; 
emeeebinspace(.4102, .4099,500); 
memes 52. te + 95.6; 

eee = .4098; 

Pammm@e = (Split split]; 

peeee = (-1.32 1.32]; 


peem~eeee- 1.32 1.32 -.6 .5]); 

memcmon plot (xline,yline, ’g-'); 

pee, VY, °C. ',aa,bdb,'’q.',cc,dd,’‘g.’,ee,ff,’g.’); 
quadrilateral 

meme, y, W.',a,b,'w.’,c;,a,'w.’,e,f,'’w.’): 

mee (  Subquad, Quadrilateral and Split of .4098’); 
ese’ about to print’) 

oho gyam 

pause;clg; 

neta off 


of oO 


oe 


oaoale s ~= 0 
rer Peers. 3 


Serinct(‘’lter=%0' ,iter) 


IRS, 


Serre as the sides are iterated 3 times. This ensures that 


falling in the subquad will not give the sequence 1100. 


Sole 


subquad 


xO(n) = 1.032" 4 eee oe 
V0 (An) =) Seno 
end 
% 1 
eonmeme = JL SON 
avin) = 1.00] 14 *aae) Je oem 
DOH 35 ie =e aa 
end 
$ 2 


CO(}) «= 20 4 ea ere 
G0: eine ese 
end 


ene e- = lee Ol 


edi k) s=71.0 =2 1 4*6{ kk) 2 
eee) = a eek: 
end 


axis (ti-1lesD Wee cee ore 


$5 = §S +1; 

(endo) criiiea | 274. lr- 

hold on; plore aera 

plot (xXx, VV ro.’ ,aa, bb, OC.” cc Cayo. eo a area eae 
quadrilateral 

plot (x0,VU0, we’ ,a0, bU, we’ FeO, dG et Ue eee 
lterated subquad 

title(’Iterated Subquad, Quadrilateral and Split of 
pause; 


joe Hie io) .25e 
Prepare for next iteration. 
X=xO0>V=V0;a=a0 > b=50 .c=cl -d-d0 =-e0 a_i. 
% ie ‘SSS 

ous UGC alue| 

end 


end; 


AG 


split 


oP oP 


oe 


4098s 


APPENDIX B: TYPICAL LEFT QUADRILATERAL RUN 
The following is a typical run of GRDCOMP6.M which 
models which points from the left quadrilateral give the 
same binary sequence as a point (W,Z) chosen at random from 
the left quadrilateral. The point (W,Z) (+) is indicated 
fee a long arrow; (u,v) points (x) are indicated with 


shorter arrows. 


Grid Field and Point Chosen on Aaractor 





Figure 55 Typical run left quad model: before iteration 
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Movement of Points in Gnd with Same Binary Sequence 
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Figure 56 Typical run left quad model: Iterate 1 


Movement of Points in Grid with Same Binary Sequence 





Figure 57 Typical run left quad model: Iterate 2 


Lae 


Movement of Poincs in Grid with Same Binary Sequence 





Figure 58 Typical run left quad model: Iterate 3 


Movement of Points in Grid with Same Binary Sequence 





Figure 59 Typical run left quad model: Iterate 4 


sls 


Movement of Points in Grid with Same Binary Sequence 





Figure 60 Typical run left quad model: Iterate 5 


- Movernent of Points in Grid with Same Binary Sequence 





Figure 61 Typical run left quad model: Iterate 6 


124 


Movement of Pounss in Grid with Same Binary Sequence 





Figure 62 Typical run left quad model: Iterate 7 


Movernent of Potnes in Grid with Same Binary Sequence 


— 


/ 





Figure 63 Typical run left quad model: Iterate 8 


dey Za) 


Movement of Poincs in Geid with Same Binary Sequence 


- 0.5 0 0.5 l 
Figure 64 Typical run left quad model: Iterate 9 


Movement of Poincs in Grid with Same Binary Sequence 





Figure 65 Typical run left quad model: Iterate 10 
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Decrease in No. of Grid Points with Same Binarv Sequence ard 4)" vs.on 
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Iteration 
Figure 66 Grid points with same sequence after n 
1terations 


Moeme 5 shows the count of points for iterates 1 through 10 
that give the same binary sequence as (W,Z). At iterate il 
(not shown) there are no grid points that give the same binary 


sequence as (W,Z). 


TABLE 5 NUMBER OF POINTS THAT GIVE THE SAME BINARY SEQUENCE 


NUMBER OF POINTS THAT GIVE THE SAME BINARY SEQUENCE 


| NO 


epee eso" |. 50 





tae 


W,Z and Grid Points with Same Binary Sequence 





Figure 67 (W,Z) and three (u,v) points prior to 
1teration : 


Movement of Pt Chosen on HA and Mauching Sequence Pt unal Diverge 





Figure 68 (W,Z) and three (u,v)’s: Iterate 1 


Movement of Pt Chosen on HA and Matching Sequence Pr unal Diverge 





“1 0.5 0 0.5 l 


Figure 69 (W,Z) and three (u,v)’s: Iterate 2 (Notice 
that the (u,v) points are essentially superimposed). 


Movement of Pt Chosen on HA and Masching Sequence Pt unal Diverge 





Figure 70 (W,Z) and three (u,v)’s: Iterate 3 (Notice 
that the (u,v) points are essentially superimposed) . 


NS, 


Movement of Pt Chosen on HLA and Matching Sequence Pt unal Diverge 





Figure 71 (W,Z) and three (u,v)‘s: Iterate 4 (Notice 
that the (u,v) points are essentially superimposed). 


Movemenc of Pt Chosen on HA and Matching Sequence Pt unal Diverge 





Figure 72 (W,Z) and three (u,v)’s: Iterate 5 (Notice 
that the (u,v) points are essentially superimposed). 


ibe) 8, 


Movement of Pt Chosen on HA and Masching Sequence Pt unal Diverge 





Figure 73 (W,Z) and three (u,v)’s: Iterate 6 


Movement of Pt Chosen on HA and Mauching Sequence Pt unal Diverge 





Figure 74 (W,Z) and three (u,v)’s: Iterate 7 


ile 


Movement of Pt Chosen on HA and Masching Sequence Pr unal Diverge 
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Figure 75 (W,Z) and three (u,v)’s: Iterate 8 


Movement of Pt Chosen on HA and Masching Sequence Pr unal Diverge 





Figure 76 (W,Z) and three (u,v)’s: Iterate 9 


eZ 


Movement of Pt Chosen on HA and Masching Sequence Pt unal Diverge 
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Figure 77 (W,Z) and three (u,v)‘'s: Iterate 10 
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APPENDIX C: 6 TO 16-TUPLE DATA 
This appendix contains the counts of n-tuples for typical 


Hénon binary sequences for n-tuples of length 6-16. 


Vv 
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| to AM. ll 
0 10 20 30 40 30 


Shifted Decimal Equivalent 





Figure 78 Incidence of 6-tuples in typical Hénon 
generated binary sequence 


Incidence 
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0 40 60 00 


Shifted Decimal Equivalent 





Figure 79 Incidence of 7-tuples in typical Hénon 
generated binary sequence 


134 
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100 
Shifted Decimal Equivalent 





Figure 80 Incidence of 8-tuples in typical Hénon 
generated binary sequence 


Incidence 


200 300 


Shifted Decimal Equivalent 





Figure 81 Incidence of 9-tuples in typical Hénon 
generated binary sequence 
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400 600 800 
Shifted Decimal Equivalent 





Figure 82 Incidence of 10-tuples in typical Hénon 
generated binary sequence : 
5 
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1000 1500 
Shifted Decimal Equivalent 





Figure 83 Incidence of 1l-tuples in typical Héenon 
generated binary sequence 


Incidence 
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1000 1500 2000 2500 3000 
Shifted Decimal Equivalent 





Figure 84 Incidence of 12-tuples in typical Hénon 
generated binary sequence 
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Shifted Decimal Equivalenc 





Figure 85 Incidence of 13-tuples in typical Hénon 
generated binary sequence 
36 
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Shifted Decimal Equivalent 





Figure 86 Incidence of 14-tuples in typical Heéenon 
generated binary sequence 


Incidence 


L3 Z 22) 
Shifted Decimal Equivalent 





Figure 87 Incidence of 15-tuples in typical Hénon 
generated binary sequence 


Incidence 
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Shifted Decimal Equivalent 





Figure 88 Incidence of 16-tuples in typical Hénon 


generated binary sequence 
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APPENDIX D: AN UNREALIZABLE 5-TUPLE 


As we saw in Chapter IV, just because a particular n-tuple 
is not found in a typical string of length 10° does nomena 
we will not tfind it in a typical string et lengt ian eer 
longer. The only way to show that an n-tuple is unrealizable 
is to prove it analytically as we have done for the sequence 
Vola On Oy. It is in this Appendix that we prove that the 
sequence {0,0,0,0,0} is also unrealizable. This specific 5- 
tuple is used in Chapter V in order.to support a conclusion as 
to why certaln sequences are unrealizable. 

The sequence (07.0, GRO, 0} as onewor the four 5-tuples that 
appeared to be unrealizable from our data in addition to those 
5-tuples that included the sequence {1,1,0,0}. These 5-tuples 
appear 1h © Ghapter 2) fake wee Consider the 5-tuple 
{0,0,0,0,0}. Suppose this binary sequence corresponds to the 


following sequence ot “x Values: "fan =e Ko hase All 


n+l 
these x values must, therefore, be at most .4098. Referring 
to equation (1.3) on page 10 the possible solutions Goma 
correspond to the region beneath the inverted parabola in 


Figure 89. Similarly, only those x,,, and x, terms less@aiiam 


Or equal to .4098 are considered (refer to Figure 90). 


Mois: 





Figure 89 Solution space of points that could give the 
sequence {0,0,0,0,0} under x,_, restriction 





Figure 90 Solution space of points that could give the 
sequence {0,0,0,0,0} under x, and x,,, restrictions 


JLB 


By containing the trapping region in a rectangle (see 
Figure 91) we see how the possible solution values can be 


further constrained. 


Minimum Rectangie to Encompass TR 





Figure 91 Trapping region in minimum area rectangle 


The minimum and maximum x values, for example, are -1.32 and 
1.33. In order to show that a string of five successive zeros 


ls not possible we consider the recurrence: 
ee = iL aa 1 aes + Bo aay (Dea 
where x, and x,., are implicit in Ene equations. Toe Optimum 


our potential of finding a value for x,,, at most x,.,=.407 sme 


require a maximum x,,, value and a minimum x,,, value. Clearly, 
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faeeeecanm be no less than -1.32. However, x Hoenoe, LimiEted 


Ne2 


only by .4098, as we shall see. Consider the equation: 


To find our desired maximum x,,, value, we must use a minimum 


wee LOr x, and a maximum value for x,.,, such that 


vee ssl nA AOS) 4 aes tag 2) 


Mies, Lurther constrains x,,, to a value no greater than 


approximately .3689 (see Figure 92). (The value of x,,, to 


elght decimal places is .36888954). 





Figure 92 Solution space of points that could give the 


sequence (0,0,0,0,0} under x,., and x,., restrictions 
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Thus, referring to equation (D.1), the maximum x,., valvieuweemno 
more than .3689 and the minimum x,., value 1S -1.32 based on 


the confines of the trapping region. It follows that 


IV 


x4 2 la) eng 2 ceo) rene asin 


x, 2 .41348871 (to eight decimal placcaim 


The value of x,,,; corresponds to a binary value Of Jae 


{0,0,0,0,0} 1s not realizable under this symbolic dynamics. 
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APPENDIX E: SUBORDINATE PROGRAMS 


This Appendix includes programs which are necessary but 
secondary to the programs in Appendix A. Two of the programs 
(HENON and HENREAL) are essential for the proper operation of 
the APPENDIX A programs. 


Fh 
G 
eo) 
9] 
ct 
bs. 
O 
eo) 
* 
aS 


F = henreal(n,x0,y0) 

NOTE: THIS PROGRAM IS REO’D TO RUN ALL GRDCOMP PROGRAMS 
This program is credited in its entirety to the author of 
reference 13;1t has been used by permission. 

program to generate n-length real sequences based on the % 
Henon 

MetSeshoe attractor. initial points fit into the 
quadrilateral of convergence described in [Hen76]. 

mhHpuUtS : 


n = length of desired sequence 
POm= Min teaal x =-value 
yO = initial y value 


Meno | teal Vector 
=mmove eal yeecror 


Hor ok KK & oO oP oP OP oO Of oP Of BO OP OP A OY OP O10 


(de = BD} 
eee v0; 
Mmeweine to check if initial points are valid 
Semewet 74 1--.1083 -1l; -3.64 1; -.1562 1]; 
Sees tll9 —-.2760 -4.6718 -.3344]’'; 
Mmmm (A~(x07v0]) > B) == 0 
disp(’‘initial point outside convergence zone’) 
meburn 
end 
S$recursive generation of points 
feomet = |:n-1; 
ae 7 ey, ed Ae (1) 2: 
Wek m= ~s "ex 1)> 
end 


peeee(xX,y,’b.'); 
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% PROBSUBQUAD .M 

% PROGRAM REQUIRES HENREAL PROGRAM 

% This program finds the dynamic probability of a palms 
falling into ehe 

% exclusion zone compared to being on the rest of the 
avEract ore 

[x,y] = henveaist too C OC ire: 


[232.7 1; =—.14 -1l; -l2 769 1]. 267 oe 
(95.6 =.29 =--3 5.251. 
= zeros (lengEn ts) 


a 

b 
l= 
Ore 


for <2 (=. lenge. 


if min((a* (x(1) 7y¥ (2) = Di == 90s Ee ie eee 

subquad 
ci Sue ee — ee 

else Ze 

end 
end 
Ae Lele a0) 
Gee evplofeicl ce 


probable ey 0 et, lengrem.) 
pause; 
subploti(2; i ds plore Veja (ee eee 
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function [xmat, ymat] = grdcomp5(sp) 

Opie, DIFF BIWN 465 IS THAT 4 CHOOSES W,Z FROM THE LEFT 
Meo THEN ENSURES 

mao, POINT IS ON THE ATTRACTOR, 5 CHOOSES FROM THE RIGHT 
UAD AT RANDOM 

meen ENSURES THAT THIS POINT IS ON THE ATTRACTOR. 

This function takes a point(W,Z)from the quadr and 
“erates the point n 

crimes in order to ensure the point is on the Henon 
tzractor(ensuring 

Moe Che nth iterate is in the left quadrant). It 
Meee ates a grid field 

based on a certain spacing, then iterates those grid field 
ints that 
match the binary string of (wnew,znew). As long as 

meee strings match, points are iterated using the Henon 
recurrence. Only that 

* \cthose) points that completely match are finally plotted 
meeording to the 
following scheme: 


of oP K_) Of IC) of OP 


of +: 


O 


& oo 'T oF H ov 


% om Seioinalecaaa £rield “points 
% w+ original (W,Z) 
S ox Original xvec,yvec point(s) which matches binary 


Seeroe Of W,Z 


format long 
% we also show how W,Z and xvec(indices),yvec(indices) 
"walk" to the Henon 
$ attractor at each iteration (with the attractor on 
screen). 
eeepc pace(-1.33,1.32,500); 
eee 4g*O + .276; 
Pemepeooace(1.32,1.245,500) ; 
meeweoe tU -4.67/18; 
g=linspace(-1.06,1.245,500) ; 
eg 5 CO -.3344; 
eee iispace(-1.06,-1.33,500); 
meer *e - 4.1119; 


Ems = . 4098; 

memes = .6314:; epi = .1e94- 

Meee -) 1314- fp2y = -.3393; 

peepee. L063 6 -1;-3.64 1; -.1562 1); 
meee - . 2/60; -4.6718; -.3344] ; 
feeeesplit:sp:1.32); 

Gree .6:50:.5); 

mem, ) = henreal (500,-1.0,-.25); 

Beeee- = [solit split]; Vilage: = i> bas pl. 2): 
Mego = length(c) ; lend = length(d) ; 
xvec = zeros(l,lenc”lend) ; ywvec = zeros(1,lenc’* lend’, 

= Oe 
Ser > = 1:lenc 
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|e pen 26 eee nic 
Ike [our marl 
%$ This next section puts everything outside of the quad to 
(2 h0e 0) 


1£ Max ((a “2 (ep) at aes) ——ane 
xvec(kK) = ep); yvec (es) —ditae 
else 
XVEC (Kiss o=- hae yvec (ae ane 
end 
end 
end 
“Vec = reshape |xvec, lene lender 
yvec = reshape(yvec,lenc*lend,1) ; 
m = £ind\(xvee == —-10 vec — lie 
xvec = xvec(m); 
YVEG = Yvecim: 
num =. vand* lenge cel. 
sWolel=p qu Rao e ehsteny ia cea VI 
Wea vec (tadex | 
Z =v Vee (index); 
Wolter % This preserves the values of W,Z 
Zulu ie 
@ige gear ao % 
ensures wattr,zattr 1s 
wattr(j+1) = 1.0 - 1.4*wattr(j)*2 + Zatti (75 6 cme 
ate rac Eel 
Za ie ee eles (ae 
if 3 > 20 & wattr(q+l) | solve sued 
end 
end 
wnew = wattr(j+1); 
Znew =a t eta) > 
wattr = wnew; 
figce CAS 16a — Ae ANG ISN. 
on Sj vec %$ preserves the values of 
KVeCCo 3 by ee 
Yy = yvec; 
axa si([—le32 1.32 =.6 a 
indices = 1:length(xx) ; S (1 *2e3 alee. ee 
= 
qvec =] lengenicc. 
qQ = indices; 
heng = s Fengt hiss) 
ho laven-: 


plot (xvec,yvec, 9g.’ ,Wattr, Zattr, Wr Fee ee, eee 


plot(o,s, “w.’,uU,V, w.',@€,f, WwW. 70,0, Weeeet oleate one 
EDZYyY - ee. 
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lod OLE; 

meme Grad Field and Point Chosen on Attractor’); 
merit ; 

pause;clg; 

while length(q) > 1; 


eeeis(([-1.32 1.32 -.6 .5]); 
P= ee emes (Si Zeiyy)) - led*xx.°2;: 
pees 3 * XK 
wO = znew +1 - 1.4*wnew.%2; 
ao. = .5*wnew: 
Reo wees Split 
Gece rime = Solat ) | x0 == split); 
else 
eee seme) < Ssolit); 
end; 
Mio lemot ml | Geemp )ya=—0. break: 
Sirse"q =  qeemp; 
end; 
ei 4 le $ counts iterates where at least 


one point matches 


qvec = [qvec,length(q)];% plot of how # with same binary 
sequence 


lenq = length(q) 
Miemees = indices(q); % decreases each time through loop 
pe= X0( qd); %$ preserves to next iterate those 
matching values of 
mae = YOlq); Sex and VO 
wnew = wQ; 
Biew = zZ0; - 
pause; 
fateeisel = ON) ; 


pater MX, VY, '9.'’,w0,z0, ’w+’, xline,yline, ‘w’); 


peueeieeres, W.',U,V,'’w.',e,f,'’w.'’,g,hn,'’w.',ftplx, fply,'o’,fp2x, 
meew © ) ; 
meee Movement of Points in Grid with Same Binary 
sequence’); 

emer t ; 

mewse; cig; hold off; 


end 
epee yl 2 3 41); axis; 
meme (1/2) .>((1:length(qvec) )-1); 


peer b> Lengten(qvec))-l,qvec,’r’, (l:length(qvec) )-1,qvec(1)* 
eee D’ ); 
title(’Decrease in # of Grid Points with Same Binary 
Beguence vs 1/2n') 

lor it: 
mause;clg; 
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% In this way, we need the index key ‘indices’ to tell to 
Which -¢ hie ep e@ timer ele 
%$ surviving iterate corresponds with respect to the original 
XVEC and yvec. 

axis ([-1.3251) 323246 oe 


lakeplie” Tejas: 
plot (xline,yline, 'w’,x,y,’bD.’,0,S,'W.’,U,V, ‘Ww.’ 7e, Lb, 7 
OW ena won, 
plot (xvec(indices) ,yvec(indices),’ogx’ ,wattr, Zattr... ee 


plot fplsjtply) Oke 2s © ae 
title(‘’wattr,zattr & Grid Pt with Same Binary Sequence’); 
preLne, 
Bavisie Clg te lene ere ee 
xnew = xvec(indices) ; 
ynew = yvec(indices); % now we show a plot of how 
wattr,zattr and the point 
% that generates the same binary sequence 
xnew, ynew 
* walk around the attractor 


xnewt = xXnew; 
ynewt = ynew; 
WatLeibl =mVichee nu: 
ie EC COM — Acie eigen 
Gere ate — eee el 
KOt. = Tae 1a enewrs, oor Vinewe: 
Vue = 3" smewe ; 
WUE = 12:0) > 14 *Weeeree 2 = aeomene > 
Z0G = 237 wore re. 
Waetru =ewte: 
ZaLert -= zm 
xMewe = oxo. 


Yoewl = 3 70in. 
axis |= lesa oe, 6. 5) 
he fa om, 


plot (x,y,’b.’,xline,yline, 'w-’,0,S,'’w.',U,V, W. ,€, =. 
ge i ah 
DIOL Capi fp ly. Oooo me ao yee 
DIOU G20e a0 Gx ote ieee ee 
title(’Movement of Pt Chosen on HA and Matching Sequence Pt 
until Diverge’) ; 
piee ier 
Dpavuseve la: 
SD LOC (0 Gl 370 See Oe ee ee 
end 
no lowenns 
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% function [xmat, ymat] = grdcomp6(sp) 

eiiay DIFF BIWN 4&6 IS THAT 4 USED WATTR, ZATTR 5 CHOOSES AT 
ANDOM FROM 

eee ee? OUADR THEN FINDS WATTR, ZATTR. 

This function takes a point(W,Z)from the quadr and 
lterates the point n 

times in order to ensure the point is on the Henon 
meeaectlor (ensuring 

Meee the Nth iterate is in the left quadrant). It 
Mieeetates a grid field 

based on a certain spacing, then iterates those grid field 
Semes that 

match the binary string of (wnew,znew). As long as 

% the strings match, points are iterated using the Henon 
meeurrence. Only that 
% (those) points that completely match are finally plotted 
mBieeorgding to the 


b- oP oP TO oP? 


kK- oP M oo 


oe 


0? 'S 


%$ following scheme: 

% om Serena lyorud field points 

% w+ original (W,Z) 

% gx QuargGimal “<vec,yveco point(s) which matches binary 
Seung Of W, Z 


moeomat long 
% We also show how W,Z and xvec(indices) ,yvec(indices) 
myc to the Henon 

peeoteereaccor at each iteration (with the attractor on 
screen). 

@emeemspace(-1.33,1.32,500); 

Sees *O + .276; 

meres lace (].32,1.245,500) ; 

weewoo*u -4.6718; 

Seeemascpace(-1.06,1.245,500); 

Mewes o5 Gg -.3344; 

@=tinspace(-1.06,-1.33,500); 


eewie = .4098; 
mee = «| 6314; OS Awe shale 
meee = —) 1314; fEp2y = -.3393; 
empe>. 4074 1; -.1083 -1; -.1562 1]; 
Meee 2. 1119; -.2760; -.3344]; 
Meee |. S23Sp:split]; 
emme-.O:SD:.5); 
Perey) = henreal (500,-1.0,-.25); 
fame = (Split split]; wabmivew= oa ale. 32 |; 
meme = length(c);: Pend. = lengihvd): 
xvec = zeros(l1,lenc*lend); Vvee = zeros(!,lenc* lend); 
K=0 > 
meee = 1:lenc 
mem © = 1 lend 
Kaeser Le 
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% This next section puts everything outside of the quad to 
(= 120 =O } 


if Max (Cal = (eno cata) |) amor no 
XVeEC(K) = GXie) Gy Veet ere 
else 
xvectkk) = -H0- 5 7vec (a=. 
end 
end 
end 
XVEC = reshape xvec,lencs lercg=ae= 
yvec = reshape(yvec,lenc*lend,1); 
m = find(xvec =="“-=l0) |S. vece-—-) ue. 
MVGC =hxvee (me 
yvec = yvec(m) ; 
num = rand*length(xvec) ; 
index “= fisn( min eee 
w = xvec(index) ; 
Z = yvec(index) ; 
Wel ea % This preserves the values of W,Z 
Zac, = 1-25 
EGieeg =) ee % 
ensures wattr,zattr 1s 
wattr(j+1) = 1.0 - 1.4*wattr(3)7*2 + zattr(4) 2 some 
aAceTa cea 
Zattr (eee oe cree tera 
Lf 7 = 20a Wwattriag!) =ssele sereak- 
end 
end 
wnew = wattr(j+1); 
ZA, =) Z elite s( eyiealee) 
wattr = wnew; 
ZAC, <= Anew 
w= KVeEC; 6 preserves the values of 
xvec & yvec 
YY = yvec; 
axis ( (=1.32 1232 =36, 25352 
ind@wees be -lengen( >, [2 eee a 
ae 
OVeGw = henge (sr 
Gu ncdices-: 
jee = leielepelau ec) 
Ne cee 
plot (xvec,yvec,’g.’,wattr, zattr, w+’ ,xline,ylane, Ven 
plot(o,s,’w.’,u,V,'’W.’,e,f, w.’,G,0, W.*, tel, tpl, © ane 
{OZ a 
(Ol Gis@t ; 


title(’Grid Field and Point Chosen Gn Ateracte. 


ESC 


jeeeint 3 
Dause;cig; 
Oole lengthi(q) > 1; 


merci (| —-iaoztl.oe—-.6 .5))3 
Pe = yy + oOnes(size(yy)) =- 1.4*xx.*2; 
Poe = =. 3 KS 
mo = zgnew + 1. - 1.4*wnew.°2Z; 
20 = .3*wnew; 
Dee OL) = ees date 
Gee toe=w Enna Oe Solis | 8x0 ==-Ssolit ).- 
else 
Gqeemow=wmrind (x0 = split): 
end; 
tf bengthiqtemp) == 0,break; 
else q = qtemp; 
end; 
=) + CL %$ counts iterates where at least 
one point matches 
qvec = [qvec,length(q)];% plot of how # with same binary 
sequence 
leng = length(q) 
micices = indices(q) ; % decreases each time through loop 
hea XO(Q); %$ preserves to next iterate those 
matching values of 
yy = yO(q); 2 “0 tand «70 
wnew = w0Q; 
Znew = zO0; 
pause; 
Meme. On; 


meee, Vy, 'g.'’,w0,z0, ' we',xline,yline, ’w’); 


MErerS, W. U,V, 'W.',e6,L, WwW.’ 79,n, w.’, fplx, iply,'o’,fp2s, 
meey, 'O'); 
mueme | Movement of Points in Grid with Same Binary 
sequence’); 

Jone alpey one 

memse- cig; hold off; 


end 
peeeemeid 2 3 4)):; axis; 
Peewee (1 /2).°>((1l:length(qvec))-1); 


meet): iength(qvec) )-l,qvec,’r’, (1: length(qvec))-l1,qvec(1)* 
meen Db’); 
title(’Decrease in # of Grid Points with Same Binary 
Sequence vs 1/2%n') 

meitnt 
Bause;clig; 
% In this way, we need the index key ‘indices’ to tell to 
which grid point the 


shape 


% surviving iterate corresponds with respect to the Giga 
xvec and yvec. 

axis{ (|—1 332) 32 2-oeeeee 
note. on- 


plot (xline,yline, ‘'w’,x,y,’b.’,0,S,'W.’,uU,V, W.’,@€,f, (eee 
oe Ve 
plot (xvec (indices) ,yvec(indices), 'gx’,wattr, Zattr,, vee 
Dict (tplx, boly 7G 4 be7 tea ones, 
title(’wattr,zattr & Grid Pt with Same Binary Sequence’); 
| eee, 
pause cla: 


ile) ols ape 2 
xnew = “Vectncdhees i. 
ynew = yvec(indices); %® now we show a plot of how 


WAaEER, ZaLeEr andere peine 
%$ that generates the same binary sequence 
xnew, ynew 

%$ walk around the attractor 


xnewt = Xnew; 
ynewt = ynew; 
Waele laaawerme ic - 
Zeb Cee ee ae 
COL eae ei ol 
xO = 10 = 14 * newt. 2a neve 
VO. 23" Snewr. 
WOE. = 190 >) dea Wwaeere. 2 + Zoe 
Z20U-:= gsreWel ee hee 
WaACELE vam 
Zatect. = Zee 
XMewL = Oia 
ynewt = yOt; 
ates — eae le ee Gy Sl le 
lava lol ieteher 


plot (x,y,’b.’,xline,yline, 'w-',0,S, ‘WwW.’ ,U,V, W.’,e, =. 
li See ie 
DLOL (ip ie ply os helen eo 
DIOE CXOE 0b, “Ox WUE ch wo 
title(’Movement of Pt Chosen on HA and Matching Sequence Pt 
until Diverge’) ; 
Digeice 
pause;clg; 
DOE Ory Cee i Cee Ole ee 
end 
igie@lioh 1128 8 


eae 


% function [xmat, ymat] = proof5(spcers,incr,spfn) 

feeg@ee ie WE USE SPCKRS TO CHOOSE A POINT AT RANDOM THEN ENSURE 
meets ON THE 

meee ee ACTOR. ET IS THEN USED AS THE COMPARISON POINT FOR ALL 
PURTHER SP VALUES 

foNbY DIFF BIWN 5&6 IS THAT 6 CHOOSES W,Z FROM THE LEFT 
@eeek THEN ENSURES 

aes POINT 21S ON THE ATTRACTOR, 5 CHOOSES FROM THE RIGHT 
QUADR. 

6 A NEW WATTR, ZATTR ARE FOUND EACH TIME SP CHANGES. 


% This function takes a point(W,Z)from the quadr and 
1terates the point n 
%$ times in order to ensure the point 1S on the Henon 
attractor(ensuring 
$ that the nth iterate is in the left quadrant). It 
initiates a grid field 
%$ based on a certain spacing, then iterates those grid field 
Bemnts that 
% match the binary string of (wnew,znew). As long as 

the strings match, points are iterated using the Henon 
ecurrence. Only that 

(those) points that completely match are finally 
considered. 


of Ko 


peers = .4098; 
Se = Spcrs; 
Zee ooo = hemo. 64 1: = 2.1562) 1); 
it] 2/60; -4.6718; -.3344]; 
MeEmsoLre sp sl.32); 
ame - OSD: .5]); 
Memes = length(c); Pee t= Lem@is mae, 
xvec = zeros(l,lend*lenc); yvec = zeros(1,lend*lenc) ; 
Kea 
memo = 1:lenc 
moe r = J:lend 
Kak] : 
Mmeamacc (a ** “fetiem-dirj) |) << b}) == 0 
aoe = Cl Dp); yvec(K)/=a(r); 
else 
pouec (ia = -)0-Yyvecit kk} = -10; 
end 
end 
end 
xvec = reshape(xvec,lenc*lend,1); 
yvec = reshape(yvec,lenc*lend,1) ; 
m = find(xvec~= -10 | yvec~= -10); 
beec = xXVeC (m) ; 
merec = yvec(m); 
mum = rand* length (xvec) ; 
mieex = fix(num) + 1; 


LSS 


wW = xveci@mdex) | 
Z.= VVC la dex). 
Wee ie aun % This preserves the values of 
W,2 
Zac Ciena Ce 
POL oe oO 
ensures wattr,zattr 1s 
wattr(j+1) = 1.0 - 1.4*wattr(q)*2 + zattr(,) > some 
avEeracecn 
Zattr( + ieee .3*watcr ee 
LES > 20k Wate et (eee eo ie ee oiece: 
end 
end 
wnew = wattr(j+1); 
ZNCW. = Zater (ie 
fOr Sp = SDinl¢ men 4 pen-- %$ MAJOR OUTER LOOP 
Sp = sp 
a-=“[ =. 1L0G3a-1) [3264 eo eae 
b = |) -32760- 4-42.67 eo ae oe 
G vee [So ist..sio he, le 
Gi = [= esp > 
lene = lenge ci) lend == lengemid)- 
xvec = zeros(1,lend*leneie. yvee == zeros lend] lua 
kee 
LOr oes seme 
fOr 2 aa aleme 
k=k+1 ’ 
if (max Cla’ * [elo eli) | ee ee 
mvecitik)) ]:¢ (Dp) 7yvec i) —ciaa 
else 
xvec(k) = —1L0;5> yvvee(k ker 
end 
end 
end 
xvec = reshape(xvec,lenc*lend,1); 
yvec = reshape(yvec,lenc*lend,1) ; 
m = find(xvec~= -10 yl vvce-— se. 
KV ECO] Seve c tm)? 
YVEC =syveeim); 
wattr = wnew; 
Zaeer = Zaew; 
oe = ee. % preserves the values of 
xvec & yvec 
yy = yvec; 
indices = 1:length(xx) ; = [| cece J 
mS be 
Gvec = Jlengempoa 
G€o=" Indices. 
lenqvec = length(xx) ; 


while length(q) 


Sy lee 
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Po = yy +wemes (Ssize(yy)) > led*xx.*2; 
vi 3 * KX; 
wO = znew + 1 - 1.4*wnew.%2; 
A=, 3 “wnew; 
ie ee = ese at 
Semen —s cme O08 = soli1t | x0 S= split); 
else 
Geemou= samen~0 =" split): 
end; 
tt Senet iweisemp) == O-break; 
else q = qtemp; 
end; 
1 =i1+ 1;%#counts iterates where 2 one point matches 
aqvec = [qvec,length(q)];% plot of how # with same binary 
sequence 
Meng = length(q) ; 
meemees = indices(qd); %$ decreases each time through loop 
Pee XO (OC) ; % preserves to next iterate those 
Waeening values of 
yy = yO(q); x0 and yO 
wnew = w0; 
aew = zO0; 


end 


° 


%$ In this way, we need the index key ‘indices’ to tell to 
Women Grid point the 


$ surviving iterate corresponds with respect to the original 
xvec and yvec. 


xnew = xvec(indices); 
ynew = yvec(indices) ; 
xnewt = xnew; 
ynewt = ynew; 
Woemrret = Wattr; 
Zee = ZAaAcCtr; 
meme = 1:1 - lL 
pee 1.0 — 1.4*xnewt.°2 + ynewt; 
Paves] .3*xnewt; 
migee= 120 - 1.4*wattrt.72 + zattrt; 
faiew=  .3o*wattrt; 


Nemert = wO0t; 
Bel@ert = 20t;> 


mmewet = x0t> 

moewt = yOt; 

end 

meee = ([nvec,n]; 

Sevee = [SPveCc, sp]; 

clear xvec; clear yvec; 

end % END TO MAJOR OUTER LOOP 
Bees - 


evot (spvec,nvec, 'w*’) 
meee’ SP VS. TTERATIONS UNTIL DIVERGENCE’ ) ; 


Les 


% function [xmat, ymat] = proof6(spcrs,incr,spfn) 
% HERE WE USE SPCRS TO FIND A POINT THEN ENSURE iT Ilse) eo 
PNM vewaS IG 

THIS POINT IS THEN USED AS A COMPARISON POINT FOR Malta 
ALUES. 

ONLY DIFF BTWN 4&6 IS THAT 4 USED WATTR, ZATTR 6 CHOOSEo a 
ANDOM FROM 

THE LEFT QUADR THEN FPINDSeweaee. 2ZA0 ce 

This function takes a point(W,Z)from the quadr and 
terates the point n 

times in order to ensure the point is on the Henon 
ti ractortensur me 

that the nth iterate is in the left quadrant). It 
niliates a G@raa. trela 

based on a certain spacing, then iterates those grid field 
Sints thar 

match the binary string of (wnew,znew). As long as 

the strings match, points are iterated using the Henon 
ecurrence. Only that 

(those) points that completely match are finally 
considered. 


o© oP "DI oO <y oO 


oP 


oe H- oP 


0 K of a 'O 


Split.= 32032; 
Sp = pers. 
a = [3.4074 1. -210e2 31 (cz 
be => [—-4ieel lo e279 60 =. 2244 
ome — el bol. 5) eae) ssye) la fe le. 
Gd. = "| 6 Spas 
bene = Vengo tiie): bend = lengencaa: 
xvec = zeros(1,lend*lenc); yvec =" zeros(1, lence" lena 
Pe 
[Or - Seeee lemme 
Reta 1S) vend 
hae. 
LEomax( (a 7 [eos cur) =o) ae ane 
XVEC(K) \= Cle) ace) —atae 
else 
xvec(k) = =05 -yveet |) Same 
end 
end 
end 
“VEC (= Meshape |(xvee, lene tena ks. 
YVGC = echape(yvec, Lene Veuce li 
m = find(xvec~= -10 | yvec~= -10); 
xVeC = xVvec(m) ; 


yvec = yvec(m) ; 


nuMe = rand“ Tenierh tayea) 
index =] fixtaum ea 
w = xvec (index) 
Z = Vvec( tadesa] 


SG 


wattr = w; % This preserves the values of W,Z 
meee r = Z; 
ao) |6U= CULL 2 AO 
%$ ensures wattr,zattr is 
Pee oe) eel Tee twattr( 7) “2 + zattr(,);% on the 
attractor 
maieer tye li = 24 wattr( 7); 
Pee ee eewartr( y+) < split, break; 
end 
end 
wnew = wattr(j+1); 
mmew = zattr(j3+1); 
mepeesoD = Spin:incr:spcrs, % MAJOR OUTER LOOP 
Sp = Sp 
ceeeeoeag74 1+ -2 1063 -1; -—.1562 1]; 
fees TIO = 2760; 980 6-33.44); 
meet i 32°Ssp-split] ; 
oe - 8 6 Sp: 5]: 
Meme = lLength(c); ene =a bend rin ci; 
Peon zeros(l,lend*lenc); yvec = zeros(1,lenc*lend) ; 
oO 
mer Oo = |:lenc 
mor r+ = 1: lend 
hehe Lee 
memmact (a “8 (clip) ;d(r) )) < b) ==.0 
Pa ee Keep yvec(k) =a lr)"; 
else 
PreGame yeci) —10; 
end 
end 
end 
xvec = reshape(xvec,lenc*lend,1); 
yvec = reshape(yvec,lenc*lend,1); 
m = find(xvec~= -10 | yvec~= -10); 
ee = xvec(m); 
yvec = yvec(m); 
wattr = wnew; 
Zaeee = ZNew; 
pee = XVeCC; %$ preserves the 
values of xvec & yvec 
VY = yvec; 
em@emees = 1:length(xx); le ren 
ah eae ee 
exec = length(xx) ; 
Co slmGiCes - 
lengvec = length(xx); 
while length(q) > 1; 
i= a7 emomes (Ssaze(yay)) =]. 4*xx.°2} 
Pee= So 
wWoP=eznew + 1 =— 1.4*wnew.*2: 
ZB =. 3 “wnew; 


eS 


Lt WO] = ic oem 


qtemp = £ind(x0 S "splat. <0) == ecw euem 
else 
qtbemp =) find (0s op maneee 
end; 
1f Jenoth (qtemp)) ==sthorea. 
else gq = qtemp; 
end; 
1=i1 + 1; % counts iterates where at 
least one point matches 
aqvec = [qvec,length(q)]; %* plot of how # with same binary 
sequence 
lenq = length(q) ; 
indices = Wheres (ee % decreases each time through loop 
<x Uta % preserves to next iterate those 
matching values of 
eve on ONben < % x0 and y0 
wnew = wQ; 
znew = 20; 
end 


° 


% In this way, we need the index key ‘indices’ to tell to 
When Gislemeeemm cine 


4 surviving iterate corresponds with respect to the original 
xvec and yvec. 


xnew = xvec(indices) ; 
ynew = yvec(indices) ; 
xnewt = xnew; 
ynewt = ynew; 
Woe bik sve el 
Zab r = Zawen> 
fOm = lal 
KOC = 1.0 = 2.4*xnewt., 2. yew ee 
V0 C= ee ewe 
wOt-= 1290 = 1.4*wattut sy 26 Zarrue. 
ZO =e Wabere. 
WaRErt) = wee 
ZabLtre = Zee 
XNewe, = ec0c; 
View: ear aac: 
end 
nvec. = Simvee a. 
spvec = [spvec,sp]; 
clear xvec;clear yvec; 
end % END TO MAJOR LOOP 
axis; 


PLlOEAS Ep VeC mMce ss wu) 
title(’SP VS; TTERATIONS UNTIL DIVERGENCE ae 


IDs 


function x = henon(n,x0,y0) 

iireoemeLregtam tS credited to the author of reference 13: it 
has been used by permission. 

am@umetion x = henon(n,x0,y0) 

program to generate n-length binary sequences based on the 
$Henon 


NP AP oO 


Shorseshoee attractor. initial points are checked %against 
the 

S$quadrilateral of convergence. (Hénon 1976) 

ee nous : 

% n = length of desired sequence 

% eUe= initial x value 

% Viw=-santeral y value 

BOULOucs : 

S Po OY ima YOCLOr 

S 

1c (| la =a Oe 

yi 1) vue 

Meett = .4098-; median x-value of henon attractor 


emeutime to check 1f initial points are valid 
pees 74 1 —. 1083 -1; -3.64 1; =-.1562 1]; 
Pee elo -,2760 -4.6718 -.3344)'; 
mumpmeer A” |xXO>yO)]) > B) == 0 
disp(’initial point outside convergence zone’) 
return 
end 


Paes) = zeros (n-1,1); $vectors are preallocated here to 


save time 
wee) = zeros(n-1,1); $1n case initial point is outside 
of zone 


$recursive generation of points 

mer 1 = l:n; 
x(1+1) 
peta |) 


convert previous point to binary 
ieee mei = SO lac 


peta) galt): 
else 
pa Os a a 
end 
end - 
Sminor housekeeping to dump the last term 
mae x (1 tn); 


ila, 


bO 


LO: 


ii 


Lee 
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